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ABSTRACT 


1, 
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This  report  documents  the  results  of  a  Small  Business  Irmovative  Research  (SBIR-Phase  10 
project  conducted  for  DARPA  focusing  on  the  use  of  all-souroe  overhead  remote  sensor  imagery  for 
monitoring  underground  nuclear  tests  and  related  activities.  This  documentation  includes:  1)  the  main 
unclassified  body  of  the  report;  2)  a  separate  ground  truth  Annex;  and  3)  a  separate  rlaiwifi^iyi  Annex. 
Autometric's  approach  was  to  investigate  the  exploiution  potential  of  the  various  sensors,  especially  the 
fusion  of  products  from  them  in  combirution  with  each  other  and  other  available  ct^teral  dau.  This 
approach  featured  empirical  analyses  of  muhisensor/multispectral  imagery  and  collateral  dau  collected 
before,  during,  and  after  an  actual  undergrouixi  nuclear  test  (named  "BEXAR").  Advanced  softcopy 
digital  image  processing  and  hardcopy  image  interpreution  techniques  were  investigated  for  the 
research.  These  included  multispectral  (Landsat,  SPOT),  hyperspectral,  and  subpixel  analyses; 
stereoscopic  and  monoscopic  information  exuaction;  multisensor  fusion  processes;  end-to-end 
exploitation  worksution  concept  developnwnt;  and  innovative  change  detection  methodologies. 
Conclusions  and  recommendations  for  furdier  R&D  and  operational  uses  were  provided  in:  1)  the 
general  areas  of  sensor  capabilities,  daubase  management,  collection  management,  and  data 

s 

processing,  exploitation,  and  fusion;  and  2)  specific  multispectral,  hyperspectral,  subpixel,  three- 
dimensional  modeling,  and  unique  unconventional  imaging  sensor  technology  areas. 
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SECnON  1 

INTBOI>UCTI<»«/RACKGR€>lJND 


Historically,  the  United  States  has  relied  on  sophisticated  seismic  monitoring  techniques  to 
detect  deep  underground  nuclear  testing  conducted  by  other  countries.  Among  the  difficulties 
encountered  were  those  of  differentiating  nuclear  testing  from  such  natural  occurrences  as  earthquakes 
or  detonations  produced  through  such  human  activities  as  mining  or  construction.  As  treaties  were 
signed  limiting  nuclear  testing  to  underground  detonations  of  150  kilotons  or  less,  deceptive  means  were 
devised  by  adversaries  to  mask,  or  reduce  the  range  of,  tmique  seismic  signals  generated  by  nuclear  tests. 
Currently,  seismic  monitoring  equipment  and  techniques  are  severely  challenged  to  identify, 
i-oambiguously,  covert  low-level  nuclear  testing.  Along  with  the  political  destabilization  of  many  of  the 
formerly  cohesive  countries  and  alliances,  has  come  the  emergence  of  several  Third  World  countries 
which  have  the  capability  to  conduct  nuclear  weapons  development 

With  political  tumult  proliferating  worldwide,  timeliness  assumes  even  greater  importance  in 
terms  of  monitoring  the  nuclear  research  ca{»bilities  of  these  emerging  nations.  The  collection  of 
seismic  nuclear  data  from  an  unsuspected  nuclear  test  site  would  indicate  that  research  is  well  along  the 
way  for  the  development  of  a  weapons  capability.  Such  a  late  warning  points  to  the  need  for  other 
means  to  provide  for  earlier  indications  of  nuclear  activities.  The  difficulty  is  that  much  of  the  initial 
research  is  conducted  in  laboratories  where  few  opportunities  exist  to  kleraify  the  extent,  duration,  or 
progress  of  offensive  nuclear  research. 

It  is  important  to  establish,  early  on,  that  nuclear  testing  has  advanced  beyond  the  laboratory 
stage  with  many  of  these  nations,  and  even  nmre  important  to  be  able  to  identify  and  monitor  further 
progress. 

During  the  summer  of  1S189,  Autometric,  Inc  responded  to  a  DARPA  Small  Business  Innovative 
Research  (SBIR)  initiative  by  proposing  a  study  entitled  ’Fusion  of  Information  from  Optical,  Thermal, 
Multispectral  Imagery  and  Geologic/Topographic  Products  to  Detect  Underground  Detonations.”  Since 
men,  the  Company  has  been  conducting  researdi  to  provide  insight  to  the  use  of  current  and  pending 
imaging  sensors  and  available  collateral  dau  to  detea  a  proliferation  of  nuclear  research  and  testing  in 
emerging  countries.  This  Final  Report,  along  with  its  Annexes,  documents  the  results  of  the  research  and 
focuses  on  means  to  supplement  seismic  collections  with  accessible  imaging  capabilities. 


SECTION  2 

PUEPOSE  Of  THE  RESEAECH 


The  purpose  of  this  SBIR  was  to  conduct  in-depth  analyses  of  current  and  pending  imagery 
collection  and  exploitation  capabilities  to  determine  their  utility  for  detecting  underground  detonations. 
This  was  accomplished  by  monitoring  a  leal-world  nuclear  event,  using  a  variety  of  imaging  sensors  and 
sute-of-the-art  data  fusion  techniques.  From  these  analyses,  a  set  of  exploiution  tools  was  developed  and 
documented,  along  with  a  workstation  concept  that  takes  into  consideration  tl^  capability  to  integrate 
selected  analytical  and  documentation  processes.  Finally,  recommendations  were  made  that  specify  the 
electromagnetic  windows  and  spectral  bands  most  useful  in  detecting  the  subde  changes  that  occur  at  the 
ground  surface  due  to  underground  detonations. 


SECTION  3 

DATA  COLLECTION  PEOGEDUEES 


Under  Phase  I  of  this  SBIR,  Autometric  investigated  several  means  for  acquiring  periodic 
imagery  over  seleaed  sites,  as  well  as  for  searching  large  geographic  areas  for  significant  changes.  Both 
classified  and  unclassified  sensor  access,  tasking,  capabilities,  and  products  were  considered  and 
procedures  were  established  to  acquire  data  on  demand.  A  classified  Annex  (2)  to  this  report  addressed 
results  from  National  assets.  This  unclassifkd  report  documented  the  tasking  of  other  collection  systems 
used.  Unclassified  satellite  and  aircraft  optical  aixl  multispearal  sensor  products,  including  Seasat  radar 
coverage,  were  investigated. 

Building  upon  the  Phase  I  methodology,  similar  Phase  II  collateral  dau  collection  processes 
were  initiated  once  a  test  site  was  selected.  The  types  of  collateral  collated  for  Phase  II  included 
historical  test  data,  topographic  and  geological  maps.  Digital  Terrain  Elevation  Data  (DTED),  historical 
and  current  weather  data,  as  well  as  historical  satellite  and  aircraft  image  coverage.  All  non-digiul  maps, 
schematics,  and  imagery  were  digitized  to  maintain  data  in  one  medium. 

Once  the  coordinates  of  the  selected  test  area  were  known,  irragery  collection  plans  were 
initiated  which  included  scheduled  satellite  sensor  tasking  as  well  as  aircraft  collections  using 
hyperspectrai  sensors  (see  Figure  3-1  and  the  classified  Atuwx  (2)).  Procedures  for  integrating  with  the 
test  directorate  and  other  test  personnel  were  established  with  the  various  Government  offices.  Details  of 
these  processes  will  be  surfaced  in  the  arulyses  of  the  products  (Section  4). 
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FIgur*  3«i  Coll«ctlon  Plan 


SECnON  4 

UNDERGROUND  TEST  PARAMETERS 


4,1  TEST  COORDINATION 

Phase  I  of  the  SBIR  established  the  feasibility  of  using  multiple  sensor  products  collected 
through  several  spectral  windows  to  detect  and  monitor  deep  underground  nuclear  tests.  Based  on  that 
analysis,  a  Phase  II  decision  was  made  to  focus  on  a  pending  U.S.  (ieep  underground  nuclear  test  at  the 
Nevada  Test  Site  (NTS).  In  November  1990,  Autometric  initiated  the  Phase  II  portion  of  the  SBIR  by 
establishing  contact  and  briefing  pertinent  organizations  who  would  be  involved  in  the  study.  Personnel 
from  the  DoE,  National  Photographic  Interpretation  Center  (NPIC),  and  the  NTS  were  apprised  of  the 
study  plans,  and  MTL  Systems,  Inc.  was  placed  under  contract  to  fly  their  hyperspectral  imaging  system 
over  a  selected  test  site  before  and  after  the  test  detonation.  With  only  a  few  weeks  to  prepare. 
Autometric,  with  excellent  support  from  both  the  DARPA  SBIR  COTR  and  DARPA  Headquarters,  met 
with  the  NTS  directorate  and  arranged  to  participate  in  the  BEXAR  test  scheduled  for  a  late  March  1991 
detonation.  The  cooperation  received  from  all  the  participants  at  the  NTS  was  outstaixling.  Without  their 
support,  much  of  this  research  would  not  have  been  possible.  Special  thanks  go  to  Mr.  Larry  Draper, 
whose  interest  and  patience  were  truly  appreciated.  Mr.  Ed  Brown  of  DARPA  was  also  of  great  assistance 
in  coordinating  with  other  organizations  supporting  this  study. 

Working  with  both  the  NTS  and  MTL  Systems,  an  operations  plan  was  produced  that  would  be 
transparent  to  the  BEXAR  operations,  yet  provide  the  data  needed  for  the  SBIR  study.  Figure  4.1-1 
contains  the  Operations  Plan  that  was  jointly  agreed  upon  and  submitted  to  the  NTS  in  February  1991. 

To  assure  other  imaging  sensor  products  would  be  available  for  use  in  the  study,  an  image 
collection  plan  was  also  produced.  Historical  imagery  of  the  site  included  Seasat  radar  imagery. 
National  High  Altitude  Program  (NHAP)  aircraft  photographic  coverage,  French  SPOT  panchromatic 
imagery,  and  Landsat  multispectral  and  thermal  coverage.  The  Landsat  system  was  tasked  to  acquire 
additional  coverage  soon  after  the  test  so  that  correlations  could  be  made  with  previous  Landsat  data. 
Arrangements  were  made  with  MTL  Systems  to  provide  video  coverage  of  the  test  simultaneous  to  their 
hyperspectral  collections.  The  NTS  was  also  asked  if  diey  would  provide  copies  of  any  imagery  that  they 
collected  during  the  test  Finally,  working  with  NPIC,  the  tasking  of  classified  resources  was  accomfdished. 
A  classified  Annex  (2)  to  this  report  provides  the  details  of  this  tasking  and  analysis. 

Since  there  was  a  need  to  establish  ground  truth  support  to  comi^ment  the  analyses  of  the 
imagery,  arrangements  were  made  for  Autometiic  persomd  to  be  oit-siie  when  the  MTL  Systems  aircraft 
was  collecting  its  data.  Among  other  ground  truth  activities,  hand-held  camera  phott^rapla  of  the  test 
site  and  surrounding  areas  were  taken.  A  Ground  Trutii  Repon  CAimex  1)  has  been  comfnled 
accordingly. 
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Oparatlons  Plan 
NTS  Collactlon 

MTL  Systems,  Inc.  will  locate  its  ffght  personnel  and  its  twin  engine  Piper  AZTEC  aircraft  (tail  No.  N40510) 
at  a  site  where  routine  coverage  of  the  selected  target  can  be  accompfshed.  The  sensor  to  be  used  it  the 
ASIS,  along  with  a  conventional  VHS  video  camera.  The  target  site  is  at  coorcfinate  N37‘’1746”- 
W116'18'46".  The  site  elevation  is  7037  feet  MSL.  The  collection  window  is  from  1  April  to  12  April, 
bradteting  the  test  and  allowing  (or  some  weather  delays.  The  first  aircraft  fight,  assuming  no  delay  in  test 
scherkile,  will  occur  on  1  or  2  April  between  1000  and  1400  hours  local  time.  This  wil  provide  coverage 
prior  to  the  test  everft.  A  minimum  of  two  near  nadir  passes  win  be  made  over  the  test  site  at  attitudes  that 
will  provide  ASIS  coverage  at  ground  resolutions  of  8  meters  and  13  meters.  The  second  flight  will 
repikate  the  first  flight  both  in  terms  of  orientation  and  altitudes  so  that  image  superpositioning  can  be 
easily  achieved.  The  second  flight  will  be  scheduled  within  24  hours  after  the  test  event,  weather 
permitting,  and  no  later  than  7  days  after  the  first  flight.  Autometric,  Inc.  wiU  provide  an  annotated  aerial 
image  of  the  test  site  to  MTL  (or  their  planning  purposes.  MTL  will  work  with  the  Government,  who  will 
specify  the  parameters  of  the  flight  in  order  to  avoid  other  sensitive  sites. 

As  soon  as  the  test  data  is  firmly  fixed.  Autometric,  Inc.,  NTS,  and  MTL  personnel  will  meet  at  the  Desert 
Rock  Nevada  Airport  on  the  date  of  the  first  fight  to  review  the  fight  plan  and  prepare  for  contingencies. 
The  fight  will  then  originate  from  the  Desert  Rock  Airport  and  return.  A  similar  approach  wiH  be  taken  for 
the  second  fight. 

Autometric  personnel  wiH  be  in  the  test  area  during  the  fl^t  to  document  ground  truth  with  a  handheld 
video  camera,  and  will  meet  later  with  MTL  and  govOTvnent  persomei  to  review  the  video  camera  results 
and  determine  the  success  of  the  fight. 

Autometric  personnel  plan  to  be  in  the  test  area  every  day  between  fights  to  become  acquairtted  with  test 
procedures.  Arrangements  win  be  made  with  Mr.  LaMffence  Draper  on  how  this  can  be  best  accomplished. 

The  MTL  aircraft  will  be  flown  Thomas  F.  Cauley,  a  certified  transport  pilot,  certification  number 
1572656.  He  will  be  acoomptmied  by  instniment  operator  Robert  L  Stevens,  SSN  121-54^573. 

The  Autometric  personnet  that  wil  be  on  site  between  15  March  and  12  April  9i  are: 

Mr.  Carrel  L  Lucas  SSN  003-12-4297 

Mr.  Mathew  Heric  SSN  252-74-3309 

Mr.  Pam  Richard  SSN  218-36-7323 

I 

The  backup  person,  should  one  ot  the  above  not  be  available,  is  Mr.  Maurice  Barry,  SSN  052-56-0708. 

Mr.  Lucas  wW  wofk  dosely  with  the  NTS  personnel  to  osure  that  the  ahcrafl  and  the  participaHon  on  the 
ground  by  Autometric  personnel  Is  transparent  to  the  test  operations. 


P^re  4.1*1 

OpwMlona  Plan 
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4.2  TEST  LOCATION  AND  DESCRIPTION 


The  BEXAR  site  was  in  a  northwestern  section  (19)  of  the  NTS.  Specifically,  BEXAR  was  atop  a 
mesa  whose  elevation  was  approximately  2134  meters  above  sea  level  (Figure  4.2-1).  The  mesa  was 
bounded  on  the  north  by  a  terrain  feature  called  Dead  Horse  Flat,  and  on  the  west  and  south  by  Siletu 
Canyon.  Geologic  collateral  data  indicated  that  the  mesa  was  comprised  of  rhyolitic  ash  flow  tuff 
underlain  by  quartz  latite  ash  flow  and  mafic/rhydite  tuff  to  a  depth  of  over  600  meters  (Figure  4.2-2). 

Vegetation  in  the  area  was  sparse,  consisting  primarily  of  stunted  cedar  and  other  evergreen 
trees/bushes.  Small  patches  of  snow  and  ice  were  observed  on  the  north  sides  of  slopes  and  clumps  of 
vegetation. 

BEXAR  ground  scars  designating  site  activity  and  road  building  were  not  evident  on  the  earliest 
historical  imagery  collected  by  the  French  SPOT  sensor  in  September  1988,  indicating  that  the  site  was 
new  and  unconuminated  by  previous  tests.  Tl»  later  SPOT  coverage  acquired  in  December  1989 
documented  initial  scarring  and  the  presence  of  drilling  equipment  at  the  site.  A  few  kilometers  away, 
three  other  abandoned  sites  appeared  to  have  similar  configurations  as  tlK  completed  BEXAR  site; 
indicating  a  series  of  related  tests  may  have  been  underway.  The  size  difference  between  the  three  site 
scars  also  indicated  BEXAR  was  the  largest  and/or  most  sofdiisticated  of  the  tests.  In  1963,  NHAP 
coverage  over  the  historical  sites  showed  no  ground  scars. 

The  BEXAR  site  was  configured  as  a  truncated  triangle  oriented  with  its  truncated  point  to  the 
north-northeast  of  the  site.  A  small  square  area  was  also  observed  attached  to  the  eastertunost  leg  of  the 
triangular  scar.  The  site  had  been  deared  of  vegeution,  bulldozed,  and  scraped  to  produce  a  level 
surface  and  appeared  to  have  been  periodically  soaked  with  water  and  tolled  to  reduce  tfaist  Two 
bermed  artificial  ponds  built  along  the  northwestemmost  edge  of  the  site  scar,  contained  a  small 
amount  of  water,  rinuned  with  ice.  A  continuous  circular  scar,  approximately  300  meters  in  radius,  with 
the  Ground  Zero  (G29  as  its  center,  was  observed  on  aerial  video  imagery  CRgure  4.2-3).  This  was  a 
defoliated  trail  containing  radiological  instrumentation  used  as  well  as  a  security  measure  to  keep 
animals  off  the  site  and  to  alert  personnel  that  they  were  approaching  a  sensitive  area. 

Oriented  in  line  along  the  southern  base  of  the  BEXAR  site  were  ten  large  metal  trailers,  a  meul 
two-story  temporary  building  and  five  other  small  trailer  configurations.  Clustered  in  the  rear  of  each 
trailer  were  collections  of  auxiliary  equipotau  such  as  air  conditioneis,  humidifiefs,  generators,  etc.  as 
well  as  personnel  vehicles  and  an  occasional  piece  of  heavy  equipment,  depending  on  the  activity  level 
of  the  ongoing  test.  Emanating  from  the  nmdimn  ends  of  the  trailets  were  huge  vidiile  cable  bundles 
which  converged  just  to  the  south  of  the  GZ,  fonning  one  huge  lot^ring  burtdle  diat  terminated  at  the 
GZ 
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4.3  TEST  CHRONOLOGY 


BEXAR  was  initially  scheduled  for  the  last  week  in  March  1991;  however,  test  preparations 
delayed  Time  Zero  (TZ)  until  3  April,  at  09:00  hours  Oocal  time).  Autometric  persoruiel  were  on  site  on 
1  April  to  be  pre-briefed  on  the  test  and  to  acquire  ground  photographs  of  the  area  prior  to  detonation. 
On  2  April,  at  the  GO/NO  GO  weather  briefing,  the  TZ  was  moved  to  09:00  on  4  April  due  to 
unsatisfactory  projeaed  wind  velocities  and  directions. 

All  activities  at  the  site  ceased  early  in  the  morning  of  the  rescheduled  test.  Wind  direction 
problems,  though,  delayed  the  TZ  an  additional  two  hours  on  4  April;  however,  a  successful  detonation 
occurred  at  11:00  hours.  The  test  produced  a  small  debris  cloud  that  moved  to  the  north  of  the  site  and 
within  2.3  minutes  it  had  dissipated  and  was  no  longer  visible  to  overhead  imaging  sensors. 

By  14:00  hours,  above-ground  activity  at  the  site  had  resumed.  Vehicular  and  personnel  activities 
were  observed  throughout  the  afternoon,  fixing  the  critical  test  timing  and  indicating  a  successful  test. 
Finally,  approximately  24  hours  after  the  test,  on  5  April.  Autometric  personnel  visited  the  site  and 
acquired  post-detonation  ground  photography.  Additional  information  concerning  site  activity  and 
conditions  shortly  before,  during,  and  after  the  test  is  documented  in  the  Ground  Truth  Report  (Annex 
1). 


4  4  HISTORICAL  TIMELINE  CHRONOLOGY 


The  United  States  has  been  conducting  contained  underground  nuclear  tests  at  the  NTS  since 
1957,  a  year  before  the  signing  of  the  Limited  Test  Ban  Treaty.  Close  to  700  welWocumented 
underground  tests  have  been  conducted  to  date,  providing  considerable  historical  textual  data  of  the 
surface  effects  produced  by  the  detorutions.  During  the  1980s,  more  emphasis  was  placed  on  acquiring 
I  imagery  of  the  tests,  and  such  data  has  shown  interesting  correlations.  The  Defense  Nuclear  Agency 

I 

I  (DNA)  has  also  been  conducting  above-ground  simulated  nuclear  tests  since  the  1970s.  Under  the  initial 

j  Phase  I  of  this  study,  these  dau  were  combined  and  used  to  establish  a  Generic  Testing  Timeline  of 

'  events  that  occur  during  test  preparation  and  conclusion. 

I  Subsequently,  the  initial  timeline  has  been  refined  and  modified  with  data  that  have  been 
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SECTION  5 

pireuminary  visual  analog  imactrt  analyses 

Both  original  analog  and  digitally-piocessed  analog  imagery  were  analyzed  in  Section  5.  Section 

7.1  provides  descriptions  on  the  specific  digital  techniques  used  prior  to  the  preliminary  visual  anlayses. 

5.1  HISTORICAL  IMAGERY 

5.1.1  Point  Positioning  Data  Base  (PPDB)  Photographs 

PPDB  photograph  numbers  27  and  28  from  control  number  77-01200,  acquired  in  calendar  year 
1977,  provided  excellent  high-altitude  stereoscopic  imagery  of  the  BEXAR  site.  (The  PPDB  is  a 
worldwide  database  but  is  rapidly  becoming  obsolete  due  to  difficulties  in  updating  the  imagery.  The 
vast  majority  of  the  database  was  produced  prior  to  1982.)  The  ground  resolution  of  the  imagery 
appeared  to  be  approximately  one  meter.  There  was  no  indication  of  ground  scarring  or  road  building 
at  the  site.  Sparse  vegetation  and  low  shrubs  or  trees  covered  the  area  at  the  top  of  the  mesa  and  a 
recent  linear  man-made  ground  scar  parallels  a  fault  line  south  of  mesa  connecting  to  the  major 
north/south  road  servicing  other  tests  to  the  noitheast  of  the  site.  The  dark-colored,  non-reflective  cover 
at  the  site  would  have  provided  a  sharp  contrast  to  any  natural  or  manmade  scarring.  The 
photogrammetric  fidelity  of  the  PPDB  imagery  allowed  for  the  topographic  mapping  of  the  area  to 
within  a  meter  in  both  relative  altitude/height  and  geographic  positioning.  This  provided  an  excellent 
database  for  future  comparative  analyses. 


5.1.2  Seasat  Radar  Coverage 

Seasat  coverage  of  the  NTS  occurred  in  1978.  Seasat  is  an  ”L*-band  synthetic  aperture  radar 
sensor  that  was  in  orbit  for  a  relatively  short  period  of  time.  Spatial  resolution  is  in  the  order  of  38 
meters  and,  under  Phase  I  of  this  research,  Autometric  detmnined  Seasat  recorded  the  existence  of 
ct^lapsed  craters  and  of  significant  road  construction.  Since  the  BEXAR  site  was  undeveloped  in  1978, 
however,  only  the  general  terrain  surface  data  was  deteoed.  If  significatu  ground  scarring  or  road 
construction  had  been  uitderway,  Seasat  would  have  likely  deteoed  it  and  would  have  confirmed 
changes  over  previous  NTS  imagery  aoqitiieil  by  a  vadeiy  of  senaors.  The  stiertgdi  of  a  Seasat-like  radar 
sensor  is  in  its  ability  to  e^KXKt  In  otUcat  stomtions;  namdy,  rqaitfless  of  time  of  day  or  weather 


12 


t 


conditions.  This  may  offer  little  in  terms  of  long-term  tests,  but  may  well  be  able  to  detect  final  test 
preparations  as  well  as  cratering  and  targe  equipment  concentrations  occurring  directly  after  a  test. 


( 


5.1.3  National  High  Altitude  Photography  (NHAP) 

NHAP  exists  only  over  the  United  States  and  some  of  its  possessions.  It  is  indicative,  however,  of 
high  spatial  resolution  photography  in  many  country  inventories  which  could  be  made  available  on 
demand.  The  latest  available  NHAP  coverage  over  the  BEXAR  area  was  flown  on  September  10,  1983. 
Two  stereoscopic  flight  lines  were  accomplished,  one  with  black  and  white  film,  the  other  with  color 
infrared  film.  Frames  123  and  124  of  black  and  white  Mission  371614-HAP-83  provided  good  stereo 
coverage  of  the  site  area  at  a  nominal  ground  resolution  of  one  meter.  The  BEXAR  site  showed  no 
scarring;  although  the  mesa  promontory  to  the  northeast  of  the  site  showed  scars  of  previous  activity 
containing  a  square  cleared  area  and  a  road  connecting  the  area  to  the  main  north/south  road.  Usually 
the  black  and  white  photography  provides  better  ground  resolutions  than  color  infrared  film  products; 
however,  in  this  case  the  color  infrared  flight  was  flown  at  a  significantly  lower  altitude  and  provided 
stereoscopic  coverage  and  ground  resolutions  of  less  than  one  meter.  As  an  example  of  the  quality  of 
the  film,  in  viewing  frames  170  and  171  of  Mission  371614-HAP-83F  color  infrared  film,  vehicles  can  be 
detected  on  the  toads,  and  the  cabs  of  large  vaits  can  be  differentiated  from  the  bodies.  Since  the  NTS  is 
arid  and  the  date  of  the  flights  was  in  September,  the  color  infrared  film  showed  minimal  ted  tonal 
information  (depicting  healthy  vegetation).  The  near  monochromatic  blue  of  the  image  provided 
limited  additional  spectral  information  beyond  that  of  the  black  and  white  film.  Original  NHAP  film 
products  provide  excellent  stereo  coverage  of  a  given  flight  line  and,  once  registered,  provide 
photogrammetric  fidelity  similar  to  PPBD  data.  Figure  5.1. 3-1  is  a  rectified  high  spatial  resolution  color 
infrared  image  of  the  BEXAR  area  prior  to  construction.  This  was  used  to  compare  with  later  coverage 
dates  to  monitor  changes  as  they  occurred. 


5.1.4  Landsat  Thematic  Mapper  (TM)  Imagery 

One  source  of  Landsat  dau  available  for  this  study  was  a  1:100,0(X)  scale  image  map  of  Pahute 
Mesa,  Nevada,  produced  in  1964  by  the  U.S.  Geological  Survey  (USGS).  The  image  map  was  compiled 
using  Bands  2,  3,  and  4  of  the  TM  dau.  The  1964  TM  covenge  site  showed  no  indication  of  activity. 
Vegeution  covering  the  top  surface  of  the  mesa  remained  undisturbed,  and  there  was  no  indication  of 
new  road  construction  to  the  site  area.  The  vegetation  appeared  to  be  lusher  dian  in  previous  NHAP 
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Figure  5. 1.3-1 
NHAP  Color  Inlrarod 


coverage,  indicating  the  seasorul  changes  recorded  through  periodic  coverages.  Because  there  were  no 
signs  of  BEXAK  activity,  it  was  decided  to  forego  acquiring  the  dau  tapes  and  to  acquire  coverage  closer 
to  test  TZ.  Figure  5.1.4-1  can  be  used  to  compare  with  later  date  TM  coverage  to  determine  the  ease  of 
detecting  test  activities  with  a  30-meter  large-area-coverage  sensor.  Examples  of  this  are  shown  in 
Section  5.2. 


5.1.5  SPOT  Panchromatic  Imagery 

The  French  SPOT  system  can  provide  stereoscopic  images  of  a  target  area  through  two 
operational  modes.  The  first,  and  most  used,  is  the  capability  of  acquiring  coverage  over  the  same  area 
from  adjacent  orbits.  The  unique  sensor  pointing  ability  thus  produces  stereo  images  that  have  a  time 
difference  of  approximately  90  minutes  and  assures  that  little  change  has  occurred  between  the  two 
images.  A  second  capability  provides  for  steieo  images  to  be  acquired  along  the  same  orbital  line,  but 
with  days  or  months  separation  in  time.  This  capability  is  used  when  unplanned  needs  for  historical 
stereo  coverage  are  encountered,  multiple  moAoscopic  coverages  are  available,  or  tasking  of 
monoscopic  data  is  requested  over  a  previoualy  covered  urget.  Such  capabilities  also  allow  for 
collecting  combinations  of  panchromatic  and  multUipectral  SPOT  imagery  over  the  urget  of  interest. 
Although  a  SPOT  image  acquires  approximaiely  tmly  10%  of  a  Landsat  scene  at  nadir,  it  is  considered 
an  efficient  collector  of  wide  areas. 

Autometric  had  SPOT  panchromatic -stereoscopic  imagery  of  the  BEXAR  area  in  its  database 
from  a  previous  unrelated  project.  The  initial  coverage  occurred  on  17  September  1988,  and  the  other 
was  collected  on  14  December  1969.  An  in-house  SPOT  stereo  mathematical  model  had  been  developed 
that  is  more  accurate  than  the  model  sold  by  the  SPOT  organization.  This  was  used  to  register  and  view 
the  hardcopy  stereo  image  on  the  Autometric  Analytical  Photogrammetric  Processing  System  (APPS- 
IV),  and  was  later  integrated  into  a  softcopy  worksution  for  additional  analyses.  In  viewing  these  data 
monoscopically,  the  15-month  time  difference  was  used  to  detect  any  changes  that  occurred  at  the 
general  BEXAR  site. 

A  cursory  comparison  of  hardcopy  images  from  each  collection  showed  that  on  the  1988 
coverage  there  was  no  activity  in  the  area  of  BEXAR;  however,  on  an  adjacent  mesa  promontory  to  the 
northeast,  a  large  cleared  and  bulldozed  area  was  observed.  The  high  reflectivity  of  the  scarred  earth 
indicated  its  recent  origin,  and  the  detection  of  large  trailers/temporary  buildings  was  a  positive 
indication  of  an  ongoing  test  The  site  was  serviced  by  the  major  north/south  road  that  joined  several 
previous  test  areas  to  the  road  network  coruiecting  the  northwest  comer  of  NTS  to  the  rest  of  the  range 
The  1969  SPOT  coverage  provided  information  that  indicated  another  test  had  been  completed  and  all 
equipment  had  been  moved  elsewhere.  The  she  scar  was  much  less  reflective,  indicating  weather,  sparse 
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Figure  5. 1.4-1 

Landsat  Thematic  Mapper  {From  USGS  Image  Map) 
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'  vegetation,  and  the  lack  of  activity  was  allowing  the  scar  to  begin  healing.  The  lack  of  an  obvious  crater 

indicated  the  test  was  conducted  deeply  underground.  The  1989  coverage  also  documented  a  new  road 
connecting  tire  southern  portion  of  the  abandoned  test  area  to  the  top  of  the  adjacent  mesa  to  the  west. 
A  large  square  cleared  area  was  observed  at  the  terminus  of  this  new  road,  and  a  drilling  apparatus  was 
detected  at  the  center  of  the  cleared  area,  llie  indication  was  the  drilling  equipment  was  moved  along 
the  new  road  either  from  the  old  site  or  across  the  old  site  from  the  north/south  road.  The  wide  radius 
turns  connecting  the  new  road  to  the  north/south  road  indicated  the  latter  to  be  the  most  plausible  case 
and  lent  credence  to  the  fact  that  the  older  site  had  been  abandoned.  Rudimentary  superpositioning  of 
the  latest  SPOT  coverage  with  aerial  coverage  of  the  BEXAR  site  documented  that  the  square  cleared 
area  was  the  position  of  the  BEXAR  GZ  and  provided  a  reliable  indication  of  the  initiation  of  BEXAR 
lest  activity.  Figure  5. 1.5-1  is  a  stereoscopic  pair  of  images  from  the  SPOT  coverage.  Using  the 
stereoscope  supplied  in  the  back  of  this  repon,  an  appreciation  of  the  terrain  and  changes  that  occurred 
between  collections  can  be  realized.  By  simply  blinking  one's  eyes  intermittently,  changes  can  be 
instantly  identified. 

5.2  CURRENT  IMAGERY 

( 

5.2.1  Current  Landsat  Thematic  Mapper  Imagery 


5  2.1.1  Prior  to  the  BEXAR  TZ 


Tasking  for  Landsat  TM  imagery  was  requested  for  a  period  shortly  before  the  TZ  of  the  BEXAR 
test.  On  7  March  1991,  the  TM  imagery  was  acquired,  aixl  tapes  were  provided  in  early  June  1991.  Figure 
5.2. 1.1-1  is  an  inuge  comprised  of  eithanced  Bands  2,  5,  and  4  from  the  above  imagery.  This  was  used 
for  a  direa  comparison  with  the  historical  1984  TM  data,  as  well  as  other  coverage  acquired  shortly  after 
the  BEXAR  TZ.  Immediately  apparent  in  the  comparison  of  this  image  with  the  historical  dau  was  the 
evidence  of  several  new  test  scars  situated  short  distances  from  the  major  road  network,  and  conneaed 
to  the  network  by  short  unsurfaced  roads  with  wide  radii  turns.  Three  of  these  scats  had  similar  shapes 
and  coiulgurations  or  were  within  a  few  kilometers  of  each  other.  The  BEXAR  scar  showed  evidence  of 
activity  by  the  highly  reflective  line  of  trailers/buildings  along  the  southern  edge  of  the  scar,  detected  as 
fluctuating  brightnesses,  and  evidence  of  water  stored  in  ponds  in  the  northwestern  cortter  of  die  site.  A 
few  kilometers  to  the  south-southwest  wu  another  highly-reflective  scar  similar  to,  but  slighdy  smaller 
than,  the  BEXAR  site.  The  high  reflectivity  of  the  soil  and  the  short  connecting  road  indicated  the  site 
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was  under  construction.  Less  than  a  kilometer  to  the  south  of  this  site  was  another  similar  site  far  less 
reflective.  This  indicated  an  earlier  test  had  occurred.  The  above  evidence  seemed  to  indicate  that  a 
continuing  series  of  similar  tests  were  being  conducted,  with  another  test  to  be  completed  within  the 
next  several  months.  Since  no  collapse  craters  were  apparent,  this  series  depicted  deep  underground 
testing. 

Another  scar  adjacent  to  and  on  a  mesa  to  the  east  of  BEXAR  had  a  different  configuration,  yet 
appeared  to  be  cormected  to  BEXAR  by  an  unsurfaced  road.  A  close  examination  of  this  road  showed  it 
crossed  the  scarred  area  and  connected  to  the  major  north/south  road  which  services  that  part  of  the 
NTS.  This  indicated  the  two  adjacent  tests  were  not  one  integrated  test,  that  the  easternmost  scar 
represented  a  completed  test,  and  the  BEXAR  test  was  serviced  by  two  roads  connecting  it  with  the 
range  road  network. 

5.2. 1.2  After  the  BEXAR  TZ 

On  8  April  1991,  the  BEXAR  area  of  the  NTS  was  again  covered  by  the  Landsat  system,  and  the 
upes  were  provided  to  Autometric  in  mid-June,  1991.  A  visual  comjjarison  was  made  between  Bands  2,  3, 
and  4  of  this  coverage  and  the  same  bands  of  the  March  1991  data.  Figure  5.2.1.2-1  depicts  the  latest  TM 
coverage.  Using  a  mirror  stereoscope,  some  indication  of  three-dimensionality  was  achieved  and  the 
detection  of  changes  was  attempted.  The  water  signature  was  still  evident,  but  the  reflective  changes  in 
the  area  of  the  trailers  were  too  subtle  to  differentiate— especially  since  the  changes  in  solar  elevation 
and  atmospheric  conditions  contaminated  the  visual  signatures.  The  conventional  use  of  TM  Bands  2,  3. 
and  4  to  conduct  visual  change  detection  was  limited  by  the  30-meter  pixel  size  of  the  TM  and  by  the 
inability  to  differentiate  changes  due  to  human  activity  from  solar  and  atmospheric  changes  occurring 
over  time.  Using  the  full  spectral  tapability  of  the  TM  arxl  adjusting  for  non-human  changes  provided 
additional  opportunities  for  documenting  subtle  spectral  changes  produced  by  human  activities. 

No  changes  could  be  observed  at  any  of  the  other  test  sites  described  from  the  March  TM 
coverage.  The  highly  reflective  scar  to  the  south  of  BEXAR  appeared  identically  to  the  March  coverage 
and  still  appeared  to  be  under  construction.  The  scar  adjacent  to  arxl  east  of  BEXAR  had  not  changed 
No  new  scarriitg  occurred  in  the  prepared  subocene. 


5.2.2 


Current  MTL  Systems  Dau 


5.2.2. 1  MTL  Coverage  Before  the  BEXAR  TZ 

On  2  April  1991,  MIX  Systems  flew  its  Airborne  Sp)ectroradiometric  Imaging  System  CASIS)  in  a 
Piper  Aztec  aircraft  over  the  BEXAR  site  at  two  different  elevations.  ASIS  is  a  spectral  radiometer  that 
collects  simultaneous  data  using  63  separate  bands  throughout  the  0.4  to  2.4  micrometer  spectral 
wiiMlow.  The  data  were  calibrated  by  MTL  Systems,  and  die  tapes  were  then  reviewed  at  Autometric.  The 
means  devised  to  review  efficiently  these  data  are  discussed  in  Section  7.4  of  this  report.  The  intent  of 
the  study  was  to  determine  which  bands  provided  the  most  reliable  information  concerning  changes 
that  might  occur  due  to  underground  detonations.  Although  hyperspectral  data  are  not  designed  for 
visual  exploiution.  Figure  5.2.2.1-1  provides  an  example  of  one  of  the  63  bands  (Band  10). 


52.2.2  MTL  Coverage  After  the  BEXAR  TZ 

The  final  coverage  of  the  BEXAR  site  by  the  ASIS  occurred  within  three  hours  after  TZ  and  because 
of  atmospheric  cloud  conditions,  was  flown  only  cxi  the  scheduled  low  altitude  flight  plaa  The  dau  were 
reviewed  at  a  digital  display  station.  Mote  detail  on  ASIS  Arulysis  is  provided  in  Section  7.4. 


5  2.2.3  Nevada  Test  Site  (NTS)  Video  Data 

Before,  during,  and  after  the  detotution,  NTS  collected  low-altitude  video  dau  of  the  BEXAR  event. 
Primarily,  these  itiuges  were  used  by  the  researchers  as  a  ground  truth  component,  and  more 
information  on  the  video  is  provided  in  Annex  1.  The  video  complimented  the  other  data  and 
provided  a  means  to  esublish  key  observable  eveiu-related  phenomenon  (Section  5.3). 
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5.3  KEY  SPECTRAL  OBSERVABLES  FOR  PENDING  DIGITAL  ANALYSES 

After  the  abovementioned  analog  data  were  reviewed,  several  spectrally-oriented  observations 
were  made  by  the  research  team.  These  served  as  the  foundation  for  the  pendiitg  digital  analyses 
conducted  by  Autometric  and  AAI,  Inc  Csubcontractcv). 

As  a  result  of  the  detonation,  surface  disruptions  were  seen  in  the  vicinity  of  BEXAR's  drill  hole. 
In  particular,  surface  material  was  thrust  into  the  atmosphere  as  a  dust  cloud,  and  after  the  event  the  dust 
drifted  to  the  northeast  away  from  the  blast  site.  Indeed,  the  surface  effects  from  the  blast  were 
concentrated  in  the  immediate  area  of  the  drill  hole,  aiKi  although  watered  down  and  trafficked,  this 
loose  spoil  material  was  relatively  unpacked  and  unstable  and  therefore  susceptiUe  to  resettling. 

It  was  apparent  from  the  data  that  the  local  surface  geology  arul  topography  influenced  the 
surface  expression  of  the  underground  blast.  The  area  of  mafor  surface  disruption  appeared  to  stop 
abruptly  at  linear  boundaries  on  the  east  and  west  side  of  the  drill  hole.  Apparently,  a  small  geologic 
surface  lineament  on  the  east  side  and  a  steep  topographic  slope  on  the  west  side  of  the  test  mesa  may 
have  restricted  and  narrowed  the  area  of  major  surface  disruption.  No  visible  surface  disruption 
occurred  on  the  outside  of  the  lineametu  or  the  steep  slope. 

A  steep  topographic  slope  occurred  off  the  tx>tthwest  side  of  BEXAR.  In  the  NTS  video,  a  linear 
cloud  approximately  300  to  500  meters  lotig  developed  above  this  steep  slope.  It  was  possible  that  the 
soil  material  on  the  steep  slope  may  have  been  loosened  and  subsided  as  a  result  of  the  blast,  and  the 
down-slope  subsidence  may  have  entrained  the  observed  dust  The  particulates  in  this  cloud  were 
slightly  darker  in  color  than  the  particulates  in  the  clouds  thrust  from  the  drill  hole  area. 

There  were  several  potential  spectral  changes  whidi  may  have  been  detected.  The  detonation 
caused  subsuntial  surface  fracturing  in  the  spoil  pile.  These  fractures  were  centimeters  across  and  less 
than  a  meter  in  length  and  were  shadowed.  Cdlectively  they  may  have  changed  die  spectral  sigtuture  of 
the  spoil  pile.  Shadowed  areas  are  illuminated  by  difluse  hemispherical  sky  radiance  rather  than  directly 
transmitted  solar  irradiance  and  thus  have  different  spectral  characteristics. 

Another  possible  change  to  the  spoil  pile  and  unvegetaied  areas  was  a  redistribution  of  particles 
due  to  vertical  displacement  Before  the  blast,  fine  particulates  in  the  soil  and  spoil  pile  were  likely 
settled  into  crado,  crevasses,  and  other  gravitational  wells.  The  larger  coarser  particulates,  which  could 
not  as  easily  be  eroded,  remained  on  the  surface.  During  die  blast,  fine  particulates  became  entrained 
by  the  verdcai  displacement  Some  were  suspended  in  the  doud  tff  particulates  and  transported  away, 
while  odiers  may  have  been  deposited  on  the  spoil  pile  surfiree.  Some  may  have  been  simply  mixed  on 
the  surface.  The  net  result  may  have  been  an  increase  in  the  fraction  of  fine  maietiai  on  the  surface. 

Another  change  may  have  been  an  erqsansion  of  die  spcMl  pile.  The  spoil  material  was 
thstrttMled  up  to  a  meter  or  hig|iK  in  depth  over  the  top  of  the  mesa  and  m  a  maimer  widdi  created  a 
radier  dOtOna  spoil  pOe  edge  with  soeae  of  die  edp  hatvhtt  aMp  stopee.  Soane  of  tfte  edges  my  have 
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been  unstable.  Most  of  the  spoil  jriks  was  traffidted  and  wateted  down  to  prevent  dust  raising.  It  was 
unlikely  the  spoil  pile  edge  was  packed  by  vehicle  trafRc  The  seismic  displacement  during  the  blast 
could  have  caused  the  sloping  edges  of  the  spoil  pile  to  "settle  out*  or  subside.  This  would  have  caused 
the  spoil  pile  to  expand. 

The  fine  particles  which  were  entrained  into  the  atmosphere  as  a  result  of  tlK  blast  traveled 
down-wind,  in  a  northerly  direction,  away  from  the  spoil  pile.  In  the  video,  the  dust  could  appeared  to 
dissipate  sto^y  as  it  moved.  It  was  suspected  that  die  fine  particles  in  diis  doud  were  deposited  onto  the 
vegetation  and  ^ound. 

It  was  likdy  that  the  surface  displacement  during  die  blast  caused  both  direct  physical  changes  to 
the  vegeution  and,  after  the  blast,  indirect  physicdogical  changes.  Potential  direct  physical  changes 
induded  uprooting  and  toppling  trees,  breaking  branches  and  lindis,  and  leaf/needle  loss.  From  a 
remote  sensing  perspective,  these  disturbances  to  the  vegetation  may  have  resulted  in  the  exposure  of 
more  forest  floor  material. 

Surface  displacements  could  have  also  caused  harm  to  the  root  systems  of  (dants  resulting  in  plant 
stress.  Typically  plant  responses  to  stress  indude  decreased  relative  turgidity,  decreased  leaf  moisture 
content  and  nutrient  supply.  Prolonged  stress  can  produce  changes  in  certain  plant  pigment 
concentrations,  such  as  chlorophyll,  carotene,  and  xanthophylls.  Most  plant  disturbances  cause  an 
increase  in  0.4  -  0.7)un  reflectance,  and  (4ant  dehydration  and  decreased  turgidity  cause  increased 
reflectance  in  the  0.5  to  2.3  region.  These  vegeul  responses  to  stress  created  by  the  blast  were  therefore 
potential  observables. 

Although  plant  physiological  re^xmses  to  the  blast  event  would  have  started  at  the  time  of  the 
blast,  it  would  have  taken  time  for  these  effects  to  become  prcmiinenL  This  time  period  could  be  as  long 
as  it  takes  for  a  given  plant  to  die.  It  was  assumed  no  radiation  from  the  blast  afSecied  die  vegetation. 

The  surface  effects  of  the  Mast  may  have  been  more  promineru  in  areas  more  susceptible  to 
seismic  displacement.  Geologic  lineaments,  for  example,  are  often  deposition  zones  of  recendy  settled 
soil  material  and  somttimes  have  sharp  tcqiografdiic  slopes.  Steep  topogra{diic  dopes  are  subjea  to 
mass  soil  movenwnts  such  as  soil  creep  and  landslide.  Massive  soil  movements  typically  bury 
vegetaUon  and  expose  bare  soO.  Dirt  roatb  and  odtet  areas  where  die  natural  soil  and  plant  material 
has  been  removed  are  sources  of  easily  erodfide  material.  Bodies  of  water  can  change  configuration, 
and  their  waters  serve  as  dqxwkoriea  for  dtt  and  ddirto. 

Each  of  diese  phenomena  were  poteittial  indicatora  of  an  underground  detomoion.  These 
observidile^,  however,  were  apecMIC  to  KXAR.  Anodim  test  site  with  different  vegeli^on  and  sod,  and 
adiete  die  tmx  ^  was  prapeied  tUferendy,  may  requife  an  eittfrely  tfifferent  set  of  observables. 
Observables  less  dependent  on  qwdBc  enyfronaamilit  cbmacieriatlGB  Csudi  m  the  presence  df  wefl- 
deSned  gibieiie  aitfiiiiee  MwHneMd^lesve  iNMmr  paeeeihl  for  untreiari  t^ppHeSbaihi.  1M  obeerv^les 
based  on  die  i^ioB  pSe  dttieiartdMtes  are  mk  tpmSHt  to  any  one  of  enittttnaieiik  tm  are 
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applicable  to  any  test  site  prepared  in  a  similar  fashion  to  BEXAR.  Covert  underground  testing  might 
rer^ire  different  observables;  paiticulariy  if  the  spoil  material  is  transported  away  from  the  drill  hc^. 


5.4  SUJ4MARIZED  UTIUTY  OF  VISUAUY-INTERPRETED  ANALOG  IMAGERY 


As  provided  in  the  previous  report  sections,  experienced  imagery  analysts  interpreted  the 
analog  dau  and  ascertained  its  overall  utility  fcx  supporting  the  detection  of  underground  testing.  The 
most  glaring  deficiency  noted  was  the  lack  of  an  up-to-date  unclassified  high  spatial  resolution  imagery 
database.  The  PPDB,  while  satisfying  the  spatial  resolution  criteria,  is  not  being  updated  and  is 
becoming  less  useful,  especially  in  urban  areas.  The  criteria  that  it  did  satisfy  was  that  of  providing  an 
exoellem  metric  database  from  which  accurate  maps  of  the  area  of  interest  could  be  produced.  The  next 
viable  sensor  product  was  the  Landsat  TM  imagery.  With  its  large  area  coverage  and  the  capability  to 
update  continuously  its  coverage  every  16  days,  it  satisfied  two  of  the  necessary  criteria  for  monitming 
global  issues.  However,  with  its  current  spatial  resolution  of  apjvoximately  30  meters,  the  products  could 
not  be  interpreted  to  produce  definitive  information  on  facility  or  test  parameters.  The  strength  of  the 
TM  was  that  it  could  record  such  activities  as  scarring  and  rmd  building  which  accompany  imderground 
test  activities.  With  continuous  global  coverage,  the  TM  should  be  used  to  monitor  areas  desigruted  as 
potential  test  sites,  and  to  document  major  changes,  such  as  road  building,  clearing  or  scarring. 
Although  hardcopy  imagery  interpreution  can  be  accomplished  if  digital  exploitation  means  ate  not 
available,  the  seven-band  spectral  dau  that  ate  collected  digitally  can  be  exploited  more  efTidently 
digiuUy,  and  the  additiorul  spectral  information  may  compensate  for  much  of  the  relatively  poor 
spatial  resolutions.  The  Seasat  radar  coverage  pointed  out  that  if  there  were  a  need  to  acciuite  dau 
during  poor  weather  or  at  night,  these  meats  could  be  used,  if  available,  to  support  such  requirements. 
However,  the  spatial  resohitiots  were  sudi  that  only  gross  changes  in  activity  would  be  delineated.  Since 
the  prepautiors  for  testing  can  be  observed  for  months,  the  need  for  immediate  coverage  assumes  a 
low  priority  (except  for  alerting  itspeccion  verification  teams  prior  to  detonatirm).  The  SPOT  system 
products  provided  better  spatial  resohitiots  than  either  the  Laixlsat  TM  or  the  Seasat  data;  however, 
each  inuge  covered  only  about  lOM  of  a  Landsat  TM  image,  and  acquires  only  three  speoial  bands  vs. 
the  seven  Landsat  TM  bands.  Its  strengths  lie  in  its  capability  to  {»ovide  stetecso^ric  coverage  of  an 
area  and  to  monitor  an  area  from  ac^oent  otMa.  The  maps  that  can  be  produced  could  provide 
accuracies  in  the  vicinity  of  10  meters  as  qsposed  to  one  meter  from  the  PPM,  However,  SPOT 
(wovides  comiiHious  global  dau  so  dm  recent  changes  can  be  recorded  as  arms  of  iniereat  evolve.  One 
other  difiicutty  enoouniered  wRh  90r  imagery  to  dm  die  ptiotidea  for  coverage  and  dtoaemination  are 
at  dm  cRsaedon  of  a  foreign  eomitry-dm  may  or  mey  not  wtoh  to  comply  widt  a  wi^iaat  or  ba  a  pan  of 
the  vet^GUfon  appeietiu.  Vtoual  interpremton  of  die  MTl,  ASUS  dMieiid  dau  in  hardcopy  wm  not  an 
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efficient  means  to  exploit  the  material.  To  view  the  large  amount  of  data,  it  was  necessary  to  use  special 
digital  tools  and  a  display  station.  These  techniques  will  be  described  in  a  later  portion  (7.4)  of  this 
report 

To  utilize  fully  the  products  investigated,  it  was  rrecessary  to  exploit  the  strengths  of  each 
product  As  an  example,  to  condua  global  test  detection  and  nx>nitoring,  the  Landsat  Ihf  would  best 
arrswer  this  requirement  Once  a  suspect  site  is  detected,  the  PPDB  data  can  provide  a  good  historical 
and  map  base  for  further  integration  of  other  data.  Should  the  PPDB  data  need  to  be  metrically 
updated,  then  the  SPOT  system  dau  can  be  used;  provided  that  the  imaging  requiremeru  would  be 
honored  by  the  country  that  sponsors  the  system  The  capability  for  classified  systems  to  suppon  such 
requirements  will  be  discussed  in  the  classified  Annex  (2)  of  this  report 

SEcn(»<r  6 

6.0  HARDCOPY  THREE-DIMENSIONAL  ANALYSES 

For  the  past  several  years,  Autometric  has  been  producing  and  installing  their  APPS-IV  for 
military/intelligenoe  community  customers.  The  APPS-IV  is  an  analytical  plotter  designed  specifically 
for  the  generation  of  digital  databases  directly  from  stereo  models  of  most  types  of  imagery.  It  has  been 
used  primarily  to  support  sensor  tasking,  urgeting,  and  support  to  the  Cruise  Missile  program  However, 
it  also  has  the  capability  for  database  generation  arxl  mapping.  Input  to  the  APPS-IV  includes  hardcopy 
transparencies  of  the  area  of  interest,  the  mathematical  model  of  the  collection  sensor,  and  propretary 
software.  Output  is  in  the  form  of  three-dimensional  viewing,  and/or  the  production  of  thematic  three- 
dimensional  databases,  including  precise  terrain  corttouring.  A  pen  plotter  produces  the  final  product 

The  decision  was  made  to  determine  whether  a  generic  analytical  plottirrg  device,  represented 
by  the  APPS-IV,  would  provide  unique  and/or  efficient  capabilities  to  support  the  detection  of 
underground  tests.  The  initial  task  was  to  analyze  the  PPDB  imagery  covering  the  BEXAR  te^  area.  Using 
an  in-house  devek^ied  math  model,  the  imagery  was  set  up  on  the  suges  of  the  APPS-IV  and  viewed  in 
’3-D*  through  its  optics.  Since  the  PPDB  was  developed  for  just  such  requirements,  an  excelleru  stereo 
image  was  perceived.  Terrain  mapping  of  the  site  could  have  been  accomplished,  but  since  this  is  a 
routine  procedure  and  the  PPDB  coUecticm  was  accomplished  long  before  BEXAR  was  initiated,  no 
formal  tiup  was  produced.  However,  a  comparison  of  the  data  was  made  whh  dm  digital  three- 
dimensional  data  acquired  from  the  90T  settsor  arxl  displayed  on  a  digital  workstation,  the  loss  in 
spatial  resbhriion  between  the  ^CTT  dtfa  and  the  PPDB  was  significant,  and  further  reaeardi  rmia  be 
conducted  on  die  miiity  of  high-Ad^  naf^i^  m  si^ipoft  die  detection  of  und^gmind  lestir^  (Since 
the  SPOT  dau  were  digitid.  kmes  hi  spatial  resolution  occurred  wtea  hardcopy  renditions  were 
generated).  Accotdln|^.  a  dedsion  was  owde  va  reutin  as  raudi  of  die  inherent  a>OT  iqwtial  resohition 
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as  possible  by  conducting  subsequent  three-dimensional  studies  on  a  digital  display  station  (see  Section 

7.5). 
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SBCTiON  7 

DIGITAL  nOCESSING  AND  ANALYSIS 

7.1  ARCHIVED  MULTISENSOR  DIGITAL  IMAGERY  PROCESSING 


Historic  multisensor  analyses  of  BEXAR  involved  the  review  of  both  analog  and  digitally- 
processed  analog  dau  of  NTS.  The  archived  imagery  was  accessed  to  provide  a  linear  history  of  the 
BEXAR  site;  the  primary  purpose  being  to  determine  what  had  occurred  at  the  spedflc  test  site  prior  to 
BEXAR.  It  was  believed  initial  changes  coukl  be  tracked  easily  since  BEXAR  was  in  a  relatively  new  part 
of  the  NTS  range  which  had  comparatively  minimal  previous  use.  This  Section  7.1  discusses  the 
processing  of  these  dau.  The  analysis  ctf  the  resultaru  processed  images  is  reported  in  Section  5. 

As  mentioned  previously,  two  sources  of  archived  aerial  photography  were  obtained.  These  were 
USGS  NHAP  arul  DMA  PPDB.  Resultant  vertical  photos  from  these  programs  ate  highly  accurate  in 
terms  of  geometric  and  spatial  fidelity,  and  because  of  this,  they  were  considered  the  most  appropriate 
source  for  creating  the  historical  database.  The  obuined  photographs  are  listed  in  Table  7.1-1. 

To  complimeru  these,  a  1964  1:100,000  three-band  Larxlsat  TM  analog  print  was  also  obtained. 

As  discussed  in  Section  5,  upon  receipt  of  these  images,  two  analyses  were  initiated:  (1)  analog 
photographic  interpreutions  and  (2)  digitally-processed  analog  interpreutions.  Each  of  the  images, 
therefore,  had  to  be  converted  to  digital  format  To  accomplish  this,  the  above-mentioned  images  were 
digitized  using  a  300  dot/inch  Howtek  Scarunaster  linear  scanner.  The  resultant  analog-to-digital  CA-to- 
D)  conversions  generated  8-bit  digital  dau  of  the  black  a;id  white  anage»  and  24-bit  (8  biu  in  each 
additive  red,  green,  and  blue  color  record)  of  the  color  images.  Subsequently,  these  dau  files  were 
transferred  to  a  PC-based  image  processing  system  (a  Compaq  386/ERDAS  7.4  system)  for  further 
review.  The  resultant  image  examples  ate  shown  as  a  composile  in  Figure  7.1-1. 

It  was  apparent  in  Figure  7.1-1  that  due  to  die  varying  spatial  resolutions  and  orientations  of  the 
original  analog  images,  interpieutkMis  would  be  cognitively  demanding.  To  address  this,  an  additional 
step  in  the  processing  coruisted  of  digitizing  a  1978  1:100,000  USGS  three-color  topographic  map  for 
purposes  of  imagery  rectification.  The  map,  in  turn,  served  as  the  geometric  truth,  and  each  archived 
image  was  rectified  to  fit  the  nap.  This  aOowed  the  georeferenced  images  to  be  omapated  more  easily 
using  image  overlay  and  nuking  techniques. 


Tri>le  7.1-1 
list  of  rtlOtM 


•  one  stereo  pair  of  NHAP  black  and  white  paper  photographs  (98123  and  98124); 

•  one  stereo  pair  of  NHAP  black  and  white  transparent  photc^raphs  (98123  and 
98124); 

•  one  stereo  pair  of  NHAP  color-infrared  transparent  photographs  (149171  and 
149172); 

•  one  PPDB  black  and  white  paper  photograph  (027); 

•  one  PPDB  black  and  white  transparent  photograph  (027); 

•  one  PPDB  black  and  white  paper  photograph  (028); 

•  one  PPDB  black  and  white  transparent  photograph  (028); 

•  one  1:100,000  USGS  three-color  map;  and 

•  one  1:100,000  three-band  Landsat  Thematic  Mapper  paper  print 
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As  a  needed  reference  in  both  this  section  and  in  Section  7.4  CHyperspectral  Analyses)  an 
overview  of  the  specific  rectification  process  is  pro^nded.  Specifically,  polynomial  equations  were  used 
to  convert  source  coordinates  to  rectified  coordinates.  This  was  accomplished  using  first-  or  second- 
order  transformations  computed  with  19  grtxmd  control  points  (GCPs).  A  least  squares  regression 
method  was  used  to  calculate  the  transformation  matrix  from  the  19  GCPs.  Since  the  aerial  photographs 
were  vertical,  or  near  vertical  a  linear  transformation  was  chosen.  As  an  example,  the  transformation 
matrix  for  a  first-order  transformation  consisted  of  six  coefiicients  (ERDAS  7.4  convention),  three  for 
each  coordinate  QC  and  Y),  as  follows: 


al  a2  a3 

bl  b2  b3 

These  were  then  used  in  a  linear  first-order  polynomial  as  follows: 

Xo  -  bl  +  b2Xi  +  b3Yi 
Yo  -  al  +  a2Xi  +  a3Yi,  i  :. 


where: 


Xj  and  Yj  -  source  coordinates; 

Xq  and  Yq  -  rectified  coordinates;  and 
al  to  3  and  bl  to  3  ■*  coefiicients  of  the  transformation  matrix. 

When  the  transformation  matrix  was  calculated,  the  inverse  of  the  transformation  matrix  was 
used  to  retransform  the  reference  coordinates  of  the  GCPs  back  to  the  source  coordinate  system.  A 
perfect  fit  for  all  GCPs,  even  using  higher  orders,  would  have  been  difficult  to  obtain;  and  therefore, 
some  error  was  acceptable.  The  accompanying  error  was  determined  as  a  root  mean  square  (RMS)  error 
and  was  the  distance  of  the  input  (source)  location  of  a  GCP  and  the  transformed  location  for  the  same 
GCP.  RMS  error  was  calculated  as  follows: 

fiWSerror  “  (XpXj)^  (YpYj)^, 

where: 

Xj  and  Yi  •  the  input  coordinates;  and 

Xf  and  Yr  -  the  retransfonned  coordinates. 
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RMSerror  expressed  as  a  distance  in  the  source  coordinates  system.  An  RMSerror  of  for 
example,  meant  that  the  reference  pixel  was  two  pixels  away  from  the  retransformed  pixel.  Accordingly, 
for  each  of  the  BEXAR  digitized  images,  the  RMSerror  set  ^  1.0  pixels,  aiKl  this  level  of  accuracy  was 
achieved  for  each  of  the  irtuge  transformation  (vocesses. 

The  final  step  in  the  rectification  process  was  to  create  the  output  file.  Since  the  grids  of  pixels 
for  the  source  images  did  not  match  the  grid  for  the  referetKe  image,  the  pixels  were  resampled,  and 
new  data  file  values  were  calculated.  These  were  calculated  using  a  nearest  neighbor  technique. 

To  determine  an  output  pixel’s  neatest  neighbor,  the  rectified  coordirutes  (Xo>  Yq)  of  ttw  pixel 
were  retransformed  to  the  source  coordinate  system.  This  was  done  using  the  inverse  of  the 
transformation  matrix.  The  pixel  closest  to  the  retransformed  coordinates  (Xn  Yf)  was  the  nearest 
neighbor,  and  the  data  file  values  for  that  pixel  became  the  pixel  values  in  the  output  image.  The  nearest 
neighbor  method  was  used  because,  unlike  bilinear  interpolation  and  cubic  convolution  processes,  it 
transferred  origirul  dau  values  without  averaging  them,  and  the  extremes  and  subtleties  of  the  data  were 
not  lost. 

The  images  originally  shown  in  Figure  7.1-1  ate  depicted  after  rectifications  in  Figure  7.1-2.  These 
images  (and  other  rectified  images  not  shown  in  Figure  7.1-2)  represented  the  A-to-D  portion  of  the 
historical  database. 

Another  element  of  the  archived  database  was  the  SPOT  imagery.  Two  SPOT  panchromatic 
scenes,  resident  at  Autometric  from  another  prior  research  projea,  were  used.  One  was  collected  on  17 
September  1988  and  the  other  on  H  December  1989. 

As  background  of  its  sensing  capabilities  for  monitoring  imderground  nuclear  tests,  a  technical 
overview  of  SPOT  is  needed.  Originally,  SPOT  was  developed  by  the  French  National  d’Etudes  Spatiales 
and  has  a  spatial  resolution  of  10  by  10m  in  the  panchromatic  collection  mode  (20  by  20m  in  the 
multispectral  mode).  The  accompanying  SPOT  system  parameters  are  given  in  Table  7.1-2 

Table  7.1-2  SPOT  Scnaor  Paraincten 


Characteristics  of 

Sensors 

Multispectral 

Mode 

Panchromatic 

Mode 

Spectral  bands  1 

0.50  -  0.59Mm 

0.51  -  0.73Mni 

2 

3 

Instrument  IFOV* 

0.61  -  0.68Mm 

0.79  -  0.89Mm 

4.13  degrees 

4.13  degrees 

Ground  Sampling 
at  nadir 

20  x  20m 

10  X  10m 

Pixels/line 

3000 

6000 

Ground  swath  wkUi 
at  nadir 

60km 

60km 

Radiometric  resolution 

8biU 

8  bio 

Image  dau  bit  rale 

25  Mb/s 

25  Mb/s 

3-1 


I 
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The  SPOT  sensing  instnunents  may  be  pointed  to  cover  ad^cent  fields  each  60km.  In  this 
configuration  the  total  swath  width  is  117km,  and  the  two  fields  overlap  by  31cm.  Since  tlw  distance 
between  adjacent  ground  tracks  at  the  equator  is  approximately  108km,  complete  eaith  coverage  can  be 
obtained  using  the  fixed  field-of-view.  It  is  also  possiUe  to  point  the  system  off-nadir,  and  sense  any 
region  of  interest  within  a  930km-wide  strip  centered  on  the  grcxind  track.  The  width  of  the  swath  actually 
observed  varies  between  60km  for  nadir  viewing  and  80km  for  extreme  off-nadir  viewing. 

If  the  sensing  devices  are  set  at  nadir-viewing  only,  the  revisit  fiecpjency  for  any  region  of  the 
world  would  be  26  days.  This  interval  may  cause  probtems  for  the  monitoring  of  undergrcxind  nuclear 
events.  Taking  into  account  the  steering  capability  of  the  instruments,  however,  during  the  26-day  pericxl 
separating  two  successive  SPOT  passes  over  a  given  point  on  the  earth,  a  point  in  (question  could  be 
observed  on  seven  differetu  passes  if  it  were  on  the  equator  and  on  11  occasions  if  at  a  latitude  of  45 
degrees. 

The  SPOT  sensors  may  also  acquire  stereoscopic  pairs  of  images,  and  this  was  the  case  for  the 
1988  and  1989  images  used  in  this  research.  The  ratio  between  the  observation  base  (disunce  between 
the  two  satellite  positions)  and  the  height  (salellile  altitude)  is  approximately  0.75  at  the  equator  and 
0.50  at  the  latitude  of  45  degrees.  Stereoaoopic  imagery  obtained  using  these  base-height  ratios  is  useful 
for  a  variety  of  photogrammeiric  appUcttions  such  as  cartographic  work  at  scales  of  1:100,000  and  map 
updating  at  1:50,000;  and  indeed,  this  capability  allowed  the  generation  of  the  digital  elevation  dau  for 
the  BEXAR  area  (discussed  in  SeokM  7.5) 

Since  the  1968  and  1969  SPOT  data  of  BEXAR  were  panchromatic  and  not  mulUspectral,  no 
sutistical  analyses  were  performed;  rather  discriminant  visual  interpretations  were  made  from  the  two 
10-meter  images.  This  included  the  produaion  of  stereograms.  Although  because  one  SPOT  scene  was 
uken  in  September  and  the  other  in  December,  the  relative  brightness  values  from  date-to-date  were 
considerably  different  due  to  variant  seasonal  energy  inputs  (and,  likely,  atmospheric  conditions). 
Accordingly,  two  sets  of  stereograms  were  generated  for  analysis:  (1)  the  first  where  both  images  were 
unenhanced;  and  (2)  the  other  where  the  September  image  was  contrast  stretched  to  3STD  and  the 
brightness  values  of  the  December  image  were  linearly  remapped  based  on  the  spectral  space  of  the 
September  image.  These  two  stereograms  are  shown  in  Figures  7.1-3  arxl  7.1-4,  respectively. 
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Figure  7.1-3 

SPOT  Stereogram  (Unonhancod) 


/:.POT  198? 

’  ■  4  ■  ■ 


SP01^:i98:5; 


Figure  7.1-4 

SPOT  Stereogram  (Enhanced) 


With  Kgard  to  the  Landsat  TM  data,  there  were  three  data  sources:  (1)  die  aforeonentioned 
three-band  analog  TM  print  which  was  digitized  and  rectified  to  die  1:100,000  USGS  map;  (2)  one  seven- 
band  digkal  scene  collected  belore  the  BEXAR  event  on  7  March  1991;  and  (3)  one  seven-band  digital 
scene  collected  after  die  event  on  8  April  1991.  While  the  three-band  image  was  discussed  earlier,  the 
other  two  seven-band  datasets  represented  the  formal  multispectral  portion  of  the  research,  and  as  such, 
they  are  discussed  in  Section  7.2 

The  MTL  ASIS  was  also  used  to  collect  data  over  the  BEXAR  site.  Flown  in  an  aircraft,  four 
scenes  were  collected  usiitg  the  ASIS:  three  before  the  event  and  one  after.  With  63  bands,  the  ASIS  data 


represented  the  hyperspectral  component  of  the 
Section  7.4 


research.  Accordingly,  these  dau  are  discussed  in 


12  MULTI^ECTRAL  PRCXXS^4G  AND  ANALYSES 


12.\  General 

For  the  imibispectial  aspect  of  die  leseardi,  two  seven-band  Landsat  IM  ^uaaeis  were  obtained 
(an  arduved  duee-lMnd  hardcopy  priia  was  used  in  Sections  5  and  7.1  but  was  not  used  here).  The 
BEXAR  event  took  (dace  on  4  April  1991;  therefore,  die  goal  was  to  determkie  if  danges  due  to  the 
BEXAR  event  could  be  detected  in  dw  (xe-  and  (MStevent  TM  imagery.  All  analyses  in  tha  Section  are 
whole-fMxel  tediniqires.  Sub()ixd  analyses,  being  a  firequeixry  domain  modeling  (trocess,  are  discussed  in 
Section  7.3. 

Before  discussing  the  digital  analysa  of  TM  imagery,  a  ladcground  of  the  Landsat  V  system  and 
its  TM  imaging  ca()abilities  for  nKMiitoring  underground  nuclear  events  is  warranted 

The  TM  sensor  on  Landsat  V  was  launched  on  1  March  1984.  The  TM  bands  were  chosen  after 
years  of  analysis  for  their  value  in  the  discrimmation  of  vegetation  type  and  vigor,  plant  and  soil 
moisture  measurements,  differentiation  of  douds  and  snow,  and  identihcation  of  hydrothermal 
alterations  in  rock  tyfjes.  TM  data  have  a  ground-projected  instantaneous-field-of-view  (IFOV)  of  30  by 
30m  (the  thermal-infrated  band  has  a  sfiatial  resolutkxi  of  120  by  120m).  System  (>arameter  infomution 
is  (xovided  in  Table  7.2-1. 


Table  7.2-1  lanrlsat  V  Thematic  Mr^per  CM)  Syatem  Parameters 


Spectral  bands 

1 

0.45-  0.52|tm 

2 

0.52-  0.60iim 

3 

0.^-  0.69Mm 

4 

0.76-  0.90|tm 

5 

1.55-  1.75Mm 

6 

10.40  -  12.50|im 

7 

Z06-  2.35Min 

IPOV 

30  x  30mfor 
bands  1-5  and  7; 

120  X  120  m  for  Band  7 

Datt  rate 

85  ySali 

Phid  quanciankw 

Sbits 

Aitilude 

705fcm 

laadHtf  V  liad  tma  Mm  of  liMlMt  l^  n,  and  Uli  ndwe  the  thaount  of 

HMdow  in  ctf  maiSf  i|i^ili>dt-«heciuiii-  dlatfcrws 


increases  the  probabilky  that  cloud  cover  wiU  be  present  over  a  ^ven  study  area.  At  an  ortnt  of  TOSkm 
(approximatdy  the  same  as  the  Space  Siuttle’s  orbiO  repetitive  coverage  raig«*g  i6  days. 

In  Section  5.3r  se^'eral  Jccy  spectral  observables  were  discussed.  Ibese  rechided  charges  lo  the 
spoil  material  at  the  immediate  BEXAR  site  (Le.,  a  series  of  small  surface  fractures,  correspoiKling 
shadowing,  paiticuhte  soil  ledishibutions,  and  paiticulaie  weathering,  and  v^etative-ielated  conditions 
from  tire  BEXAR  site  and  surrounding  areas.  The  question  at  this  point  of  the  researdi  was:  once  the 
atmospheric  and  natural  changes  re^dertt  in  the  TM  datasets  were  removed,  could  the  event-related 
changes  be  detected  using  the  dale-ioKiate  imagery  sets.’  Thus,  before  inve^igating  diese  changes,  the 
changes  unrelated  to  the  evera  had  to  be  analyzed  and  removed  as  much  as  possible. 

Most  non-event  dianges  in  the  imagery  were  generated  by  differitrg  energy  input  quantities  and 
associated  vegetative  cydical  states.  Indeed,  (hiring  the  firra  TM  scene  coiOeaed  on  7  March  1991,  the 
solar  elevation  an^  was  37  degrees,  and  the  solar  azimuth  angle  was  136  degrees.  Con^iaratively,  during 
the  second  collection  on  8  April,  the  solar  elevation  an^  was  48  degrees  and  the  s(^  azimuth  was  129 
degrees.  There  are  a  variety  of  meduxis  for  removing  the  impacts  of  variant  energy  and  cyclical  sutes; 
and  to  enable  accurate  date-todate  change  derections,  some  appropriate  mediod  mist  be  considered. 
Two  methcxls  were  used  to  achieve  this;  one  by  Autometric,  arxl  die  odier  by  the  contraa  subcontraaor 
Applied  Analysis,  Inc  CAAI).  Because  the  validities  of  the  pending  analytical  sections  are  dependent  on 
the  meduxis  of  atmospheric  removal  employed,  diese  two  meduxis  and  their  accompanying  rationales 
are  presented  separately  in  the  following  sectkxis. 


7.2.2  Automeiric  Atmospheric  Corrections 


Because  the  solar  elevaticm  angle  of  the  post-event  TM  image  was  11  degrees  higher  than  the 
pre-event,  more  input  energy  was  reaching  the  BEXAR  test  site  on  the  second  c(dlection  day.  With  this 
increased  energy,  the  BEXAR  sRe  and  surrounding  areas  reflected  this  higher  level  back  to  the  sensor, 


and  accordingly,  BEXAR  had  a  higher  relative  reflectivity  during  the  second  collection.  To  test  for 
spectral  changes  related  to  die  event,  then,  the  energy  input  levels  tor  the  post-evmit  ct^ections  had  to 
be  modeled  and  scaled  to  matdi  more  atxurately  the  energy  input  levels  of  the  pre-event  collections. 

The  prinrary  difficulty  in  determining  valid  and  reliable  values  of  atmo^riimic  paratnete*s  such 
as  inacflarKe  and  spectral  ratflanoe  Is  in  acoouraing  for  the  decDomagnetic  energy  scmtmed  more  than 
once  as  it  passes  through  the  itmosfdiefe.  Generally,  though,  electromagnedc  etimgy  reaching  the 
ground  is  comprised  of  two  coinfnMnn.  The  first  is  dtocct  dectronragnetic  eneqy  (Emm  die  surO  whuii 
passes  ^ough  the  aiasosphe«e  whhout  being  rdiaoibed  or  scMsaeed  The  seoMid  4s  Mmospheric,  or 


i  tt  least  The  spesBd 


given  surfime  of  die  ground  is  dm  sum  of  ffie  hcadianoes  of  the  dbect  etedronraipeiic  enmgy  and  of 


diffuse  electromagnetic  energy.  The  spectral  irradianoe  of  diiea  energy  can  be  estimated  by  the  zenith 
an^es  of  the  sun,  assocated  optical  path  lengthi,  and  monpchrormuic  solar  constants.  The  leOectance 
of  a  sur&ce,  dimi,  is  die  nttib  erf  die  upward  to  downward  fluxes  of  radiaru  energy  at  the  sur&ce.  Impacts 
from  diese  energies  combine  with  spectral  compositions  of  terrestrial  phenomena  to  define  the  image 
histogram. 

The  first  step  in  removing  the  ^fects  o(  die  two-date  imagery  involved  a  review  of  the  respective 
itnagp  )stari.«sfirs.  To  do  this,  subscenes  centered  over  the  BEXAR  test  site  were  extracted  for  each  date, 
and  initial  sutistics  were  generated.  These  are  provided  in  Table  7.2. 1-1.  As  shown,  there  were 
considerable  differences  found  in  the  date^o-date  aimparisons.  From  the  sutistics,  it  was  apparent 
some  dianges  in  the  imagery  were  associated  widi  a  comUned  effect  of  increased  energy  inputs  and, 
likely,  vegeutive  cyclical  states.  For  instance,  the  difference  among  means  for  Band  5  (30  suggested 
there  were  increased  turgidity  levels  in  the  area  plants.  Indeed,  this  was  probably  the  result  of  plants 
absorbing  water  from  recent  snow  melts  (from  a  mid-March  1991  stkjw  storm).  The  increases  in  Bands 
2,  3,  and  4,  however,  suggested  the  plants,  although  possibly  having  absorbed  the  water,  have  not 
increased  their  chlorophyll  absorption  levels  appreciatively.  This  could  have  been  due  to  the  following: 
1)  the  vegeution  of  the  area  was  evergreen,  and  its  natural  cyclical  chlorophyll  absorption  levels  may 
not  have  changed  conskferably  within  the  time-frame  of  the  two  itiuging  dates;  2)  the  vegeution,  having 
adopted  to  the  harsh  desert  environment,  mainuined  a  relatively  consunt  level  of  chlorophyll 
absorption  throughout  the  period;  and/or  3)  the  vegeution  had  not  utilized  the  absorbed  water 
completely  because  night-time  temperatures  in  early  April  remained  too  cold,  and  this  prevented  the 
vegeution  from  entering  its  rutural  seasotul  growth  patterns. 

In  comparison  to  apparent  vegeutive  dunges,  however,  Band  7  suggested  that  a  substantial 
amount  of  the  spectral  reflectance  increase  was  nonetheless  due  to  the  faa  that  mote  etwrgy  was 
teaching  the  BEXAR  test  site. 

To  remove  the  non-event  changes  from  the  post-event  TM  imagery  a  series  of  histogram 
remapping  modeto  were  employed.  Specifically,  die  histograms  from  the  post-event  dau  were  expanded 
so  the  lowest  occurring  value  became  0,  and  then  the  histogram  was  linearly  scaled  so  the  maximum 
value  became  255.  Thus,  this  htsiogtam  expansion  process  took  advantage  of  the  full  8-t^  dau  range. 
Associated  diift  and  scale  fectoci  were  then  cakailaied  based  on  the  cmnposkkm  of  the  pre-event  and 
expanded  post-event  dau.  For  instanoe,  given  a  calculated  shift  fector  of  seven  and  a  calculated  scale 
factor  of  0.68,  a  post-event  pbasTs  bt^htness  vahie  (XD  would  be  re-mapped  using  the  eqm^  XI'  -  XI 
•  0.68  +  7.  The  process  resKwes  cernddenb^'  the  direa  impacts  of  inemued  ene«y  taputs  as  well  as 
radiance  and  shadowii^  impacts.  Original  pre-event,  odgitul  post-event,  and  corrected  post-event 
Mstistics  are  |»ovided  in  Tkble  7.2.1-2,  and  an  associated  unetrfiaiioed  image  aiduet  fbr  eadi  of  the 
correqwndtng  bands  is  diosvn  in  Flgaie  7.2.M. 
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7.2.3  AAI  Atmospheric  Corrections 

As  a  subcontractor,  AAI  was  part  of  the  research  team.  Their  primary  goal  was  to  implement 
multispectral  subpixel  processing  techniques  for  detecting  event-related  changes.  This  process  is 
discussed  in  more  detail  in  Section  7.3.  The  ultimate  subpixel  research,  however,  was  built  upon 
multispectral  and  hyperspectral  portions  of  the  research,  and  components  of  these  were  performed 
both  in  conjunction  with  and  independeruly  from  the  Autometric  laboratory  processings.  Acccvdingly, 
AAI  has  developed  an  atmospheric  correction  modeling  technique  within  the  domain  of  their  AASAP*. 
As  the  foundation  for  subsequent  complementary  AAI  subpixel,  multispectral,  and  hyperspectral 
sections,  their  method  of  cmrecting  for  the  daie-to-date  collection  variations  is  presented  here. 

As  ipdirated  previously,  several  multispectral  dianges  were  produced  by  the  changes  in  the 
solar  elevation  and  aaimuth  angles  between  the  pre-  and  post-event  scenes.  Again,  the  TM  scenes  were 
used  for  this  study  were  measured  on  7  Mardi  1991  and  8  April  1991,  and  the  solar  elevation  angle 
increased  from  37  degrees  (pre-event)  to  48  degrees  (post-evenO-  Similarly,  the  solar  azimuth  angle 
changed  from  136  degrees  to  129  degrees. 

Among  the  effects  that  resulted  from  change  in  the  solar  elevation  angle  was  a  change  in  the 
relative  contribution  of  radiance  from  the  terrain  and  atmosphere.  Irxleed,  the  incident  irradiance  was 
less  attenuated  by  path  absorption  and  scattering  in  the  8  April  scerre  because  solar  radiance  traversed  a 
shorter  path  through  the  atmosphere.  Reflected  radiances  from  the  terrain  and  the  atmosf^ric  colunm 
above  the  terrain  were  consequently  larger  in  the  post-event  scene.  The  contribution  from  the  terrain 
increased  more  than  that  from  the  atmosphere,  however,  due  primarily  to  two  effects.  First,  the  fraction 
of  incident  solar  irradiance  scattered  from  the  terrain  in  the  direction  of  the  sensor  was  controlfod  in 
part  by  the  bidirectional  reflectance  distributirm  functions  of  the  surface  materials,  which  had  different 
sun  angle  dependence  than  that  scattered  by  the  atmosi^te.  For  diflusely  scattering  materials  lying 
horizontally  on  the  surface  (aitd  for  a  nadir  viewing  sensor),  the  fraction  of  radiance  scattered  from  the 
terrain  toward  the  sensor  increased  in  proptMtkm  to  the  change  in  the  cosine  of  the  s<^  zenith  angle. 
Second,  there  was  an  increase  hi  the  fraction  of  directly  illuminated  (vs.  shadowed)  fortain  widi  increase 
in  sun  an^.  The  laoer  effea  had  a  more  complex  angular  dependence,  but  in  general,  was  a  smaller 
effea  than  the  former.  The  estimated  combined  effect  produced  about  a  35%  Increase  in  terrain 
radiance  frmn  7  March  to  8  April,  based  on  radiation  transfer  modd  cakufetions  using  the  AASAP.  In 
contrast,  the  contribution  from  the  atmospheric  radiance  increased  by  a  signlfkandy  smaller  amount 
due  to  its  strongly  forward-scaneiing  character.  The  increue  in  scanering  was  wavelengih  dqiendent,  but 
it  wis  in  general  less  than  about  14%,  based  on  LOVniAN  7  model  caloilatkms.  The  spectral 
diaracieristfcs  of  the  terrain  and  the  atmosiriwric  radianoe  contrftMtioos  differed  stgiifficandy,  and  the 

*  AASAP  fe  proprietary  to  AAI,  Inc 
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fact  that  the  relative  contributions  of  the  two  sources  changed  from  7  March  to  8  April  implied  that 
some  spectral  change  between  the  two  scenes  was  attributed  to  evettt-utuelated  eGfi^ 

A  second  event-uruelated  spectral  change  was  attributed  to  changes  in  spectral  duracteristics  of 
the  terrain  radiautce,  due  to  chauiges  in  the  extent  of  shadowing  from  7  March  to  8  April.  The  sunlit 
fraction  of  the  surface  was  illuminated  by  both  directly  transmined  solar  irradiance  and  diffuse 
hemispherical  sky  radiance,  while  the  shadowed  fraction  was  illuminated  by  hemisphericai  sky  radiautce 
alone.  The  sunlit  fraction  was  thus  not  only  brighter,  but  it  had  different  spectral  diauacteristics  thaui  the 
shadowed  fraction  due  to  the  differences  in  spectral  characteristics  of  the  two  illumination  scenarios. 
The  increase  in  solar  elevation  angle  from  7  March  to  8  April  decreased  the  fraction  of  the  visible 
surface  that  was  in  shadow,  yielding  a  change  in  the  effect  source  illumination  spectral  characteristics. 
This  would  have  produced  a  change  in  terrain  radiance  spectral  characteristics,  even  if  the  reflectance 
spectra  of  the  terrain  material  had  stayed  the  same. 

Both  of  these  effects,  namely  the  change  in  relative  contribution  of  terrain  and  atmospheric 
radiance  and  the  change  in  extent  of  shadowing,  were  modeled  by  AASAP  using  radiative  uansfer 
approaches.  Their  contributions  to  the  observed  change  were  consequently  identified  and  potentially 
isolated  from  other  changes  in  the  scene. 

Other  natural  changes  unrelated  to  the  event  were  not  so  easily  modeled.  These  included 
changes  in  the  density  and  growth  stage  of  vegetation,  changes  in  sky  radiance  characteristics  due  to 
differences  in  haze  and  cloud  conditions,  changed  produced  by  rain  and  other  weather-related  events 
that  may  have  occuned  between  scenes,  and  changes  produced  by  human  activity,  among  others. 
AASAP  included  a  series  of  atmospheric  and  sun  angle  corrections  as  enabling  components  of  its 
scene-to-scene  robustness  feature.  Using  header  and  relevant  elevation  formation,  corrections  were 
generated  that  attempted  to  compensate  for  differences  between  scenes  of:  1)  atmospherically  scattered 
solar  radiance-,  2)  atmospheric  path  attenuation  of  directly  transmiaed  solar  irradiance;  3)  diffusely- 
scattered  hemispherical  sky  radiance;  and  4)  fractions  of  directly  illumiiuted  vs.  shadowed  surface. 
AASAP  attempted  to  bring  scenes  into  environmental  equivalence  to  a  reference  scene,  which  was  used 
to  derive  the  signature  of  a  noaterial  of  interest.  For  this  research,  the  pre-event  (7  March  1991)  scene 
was  considered  as  a  reference  scene  and  the  post-event  (8  April  1991)  was  brought  into  environmental 
equivalence  with  the  reference  scene  using  correction  feaors  generated  by  AASAP. 

The  first  environmental  correction  applied  was  subtraction  of  AASAP-generated  atmospheric 
radiance  components  from  both  scenes.  The  corrections  are  listed  in  Table  7.2.2-1.  They  were  applied 
by  subtracting  the  listed  image  plane  data  number  from  each  pixel  in  the  corresponding  scene  image 
plane. 

A  second  etivironmetual  correction  factor  was  applied  to  the  post-event  scene  atxi  attempted  to 
compensate  for  diHierences  in  atmospheric  path  attenuation,  diffusely  scattered  hemispherical  sky 
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ndianoe,  and  shadowing.  The  AASAP-generated  correction  is  shown  in  Table  112.-2  and  it  was  applied 
ffluhiplkatively  to  die  image  plane  pixel  data  nundieis  after  subtraction  of  die  Tat^  lH-\  cmiection. 

Apiriication  of  the  corrections  in  TaMes  7.2.2-1  and  7.2.2-2  btou^  the  8  A{»il  1991  scene  into 
appraximate  environmental  equivalence  to  the  7  March  1991  scene.  Local  scale  differences  caused  by 
terrain  slope,  macroscale  shadowing,  and  topogtaphk:  elevadon  were  not  compensated  for  by  these 
corrections.  OifFeimces  in  atmospheric  haae  profiles  were  also  not  compensated.  To  oimpensate 
partially  for  these  effects,  pixels  were  selected  from  the  pre>  and  post-event  scenes  judged  to  have 
equivalent  terrain  slope,  minimized  nacroscale  shadowing,  and  topc^raphic  elevation.  Selected  pixels 
also  avoided  vegetated  areas  to  minimize  differences  caused  by  variations  in  vegeution  density,  growth 
cycles,  and  species  diversity.  Reference  pixels  were  selected  from  areas  Judged  not  likely  to  be  changed 
by  the  detonation  event  to  develop  a  compensation  for  the  differences  in  atmosfrfieric  haze  profrles, 
and  ocher  inaccurately  compensated  envirotunentai  faaors. 


Tabi*  7.2.2-1 
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Thematic  Mapper 
Image  Plane 

7M«t^1991 

CofrBcaonjtDM) 

8  April  1991 
COTTBCliOOilM 

1 

47 

53 

2 

15 

17 

3 

11 

12 

4 

5 

5 

5 

1 

1 

7 

1 

1 

Tabla  7.2.2-2 
Sun  Angla  Corractlona 
(umtlaaa  MuNlpHeatlon  Factor) 
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In  TaUe  7.2.2-3  aie  pmemed  two  sets  of  four  pbids  each,  sdtecied  6om  identical  areas  (within 
the  accutacy  of  scene  re^tiration)  in  die  pre-  and  post-event  scenes.  ‘Iheae  were  jutted  to  be  far 
enough  away  firtMn  the  detonation  ^  so  as  to  be  Uhdy  unafTetxed  by  the  detonation.  Ihe  two  sets  of 
pixels  were  corrected  using  the  foctois  in  Tables  7.2.2'!  and  7.2.2-2.  The  two  sets  corrected  pixds  were 
then  averaged  to  create  mean  reference  pixels  for  the  pre-  and  post-event  scenes.  The  pre-event  scette 
reference  pixel  was  then  divided  by  the  post-event  scene  reference  pixel  Omage  {^ane  dau  number  by 
image  plane  data  number)  to  derive  the  haae  correction  multiplier  shown  in  Table  7.2.2-4.  The 
cmiection  suggested  there  was  a  hitter  atmo^heric  haze  content  in  die  8  April  1991  post-event  scene; 
although  it  is  uncertain  where  in  the  profile  the  hitter  haze  content  occurred.  The  haze  correction 
faaor  was  applied  multiplicatively  to  the  post-event  scene. 

After  all  corrections  were  applied  to  die  post-event  data,  subscene  image  statistics  were 
generated  for  the  corrected  six  TM  bands.  These  ate  shown  in  Table  7.2.2-5  together  with  original  pre¬ 
event  subscene  statistics. 


Table  7.2.2-3 
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Table  7.2.2-4 

Haae  Correction 

Themailc  Mapper 
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Tabi»7.2^ 

kmg*  Statistics  for  ThMiwlie  ItapfMr  (TM)  SubtoMM 
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lO 

1 

CD 

Band  7 

MIN 

21 

14 

26 

34 

58 

31 

MAX 

97 

62 

109 

100 

184 

110 

MEAN 

37.25 

25.37 

45.94 

51.72 

97.38 

57.06 

STD 

10.59 

7.41 

13^4 

9.72 

19.03 

11.92 

7.2.4  Autometric  Multispectial  Analyses 

Radiometric  resection  defines  the  sensitivity  of  a  detector  to  differences  in  signal  strength  as  R 
records  the  radiant  flux  reflected  Cot  en^ttsd)  fion  the  terrain.  Radiant  flux  recorded  by  dre  TM  sensor 
records  data  in  8  (0  to  255).  Using  this  rai^e,  die  overwhdming  majorRy  of  digRai  image  analyses 

are  dependent  primarily  oti  the  resulting  muhispectral  tones/colcus  of  obiects.  Usually,  these  are 
accomplished  using  fundamental  statfetical  pattern  recognition  techniques,  (jenerally,  automated 
devices  cannot  reliably  process  elements  such  as  ^ape,  texture,  and  pattern,  let  alone  the  more 
complex  elements  of  site  and  association;  although  advances  are  being  made. 

Given  a  perspective  of  digRai  image  processing,  the  goal  of  the  muhispectral  portion  of  the 
research  was  to  utilize  the  8-bR  radiometric  resedution  of  the  pre-event  and  corrected  post-event 'm 
datasets  and  generate  models  to  detect  event-related  changes. 

The  selection  of  appropriate  change  detection  models  was  based  on  several  factors.  First,  the 
physical  characteristics  of  the  study  area  and  type  of  potential  change  had  to  be  considered.  For 
instance,  in  some  other  unrelated  area  where  change  has  been  dramatic  and  multitopical,  image 
differencing  or  ratioing  techniques  would  reveal  and  help  quaittify  subjea  classifications  (e.g.,  a  topical 
change  from  water  to  land).  Any  change  at  BEXAR,  however,  would  not  have  been  such  a  contrast  or 
multitopical  change,  and  it  would  have  been  difficult  to  differentiate  the  significance  of  such  slight 
changes  in  the  arithmetic  images.  Second,  the  data  to  be  compared  must  be  registered  to  a  defined 
registration  precision  range.  Indeed,  because  of  Landsafs  predictable  orbR,  the  two  IM  images  were 
registered  to  less  than  1.0  pixel  by  determining  X  and  Y  upe  read  parameters;  in  turn,  no  directional 
and/or  scale  rectifications  were  needed.  Finally,  knowledge  of  the  amouiu  of  diaitge  must  be  considered 
on  a  unitopical  basis  (e.g.,  did  the  chlorophyll  absorption  level  of  a  given  tree  change  from  one  date  to 
another?).  If  a  unitopical  change  had  been  great  enough,  a  vector  change  analysis  in  multidimensional 
spectral  space  would  have  been  appropriate.  The  change  af  BEXAR,  however,  was  likely  not  great  enough 
to  be  detected  using  this  method. 

Based  on  these  considerations,  the  most  appropriate  modeb  seemed  to  be  those  based  on 
investigations  of  spectral  variance.  Since  these  variations  from  date-to-date  would  be  slight,  statistical 
confidence  levels  were  needed  to  support  any  modeling.  Three  techniques  were  used  to  investigate 
variability  and  spectral  change  of  BEXAR:  1)  a  series  of  comparisons  of  band  pair  means  (Le.,  BatKl  1 
pre-  vs.  Band  1  corrected  post-event.  Band  2  pee-  vs.  Band  2  corrected  post-event,  etc.)  to  determine  if 
differences  between  bands  were  statistically  significant;  2)  a  one-way  analy^  of  variance  CANOVA)  to 
evaluate  differences  between  afl  pib-  and  oonected  post-evem  TM  bsrxis;  and  3)  a  posteriori  Tukey 
Honestly  Significant  Difference  CHSD)  test  to  eaiabiith  a  iflffHenoe  matrix  for  all  barxl  means. 

On  a  band-by-band  basis  (exchitfing  thermal)  spectral  change  observables  cm  be  estimated. 
For  instance.  Band  1  C0.4S  -  0.52  |tm  Cbltie))  provided  inforeuKton  oti  water,  soil,  and  v^etation 


SO 


itf^  WiMieaiiai  onnm  ^  Barit  of  ^  ddoro^i^  a^Mrptkxt  for  heaUiy 
^nigiStM.^m-  i  ~  t^ioB  bearoen  Uue  aid  red  dikMoptiyll 

absorption  »fld  comg^oadtet  eo  tit  lietidif  n^pnoion.  Band  2  ;d8o  provided 

i^i»oiatk»  lAiont  water.  Band  $  -■  0.(^ttift  On®  aeatiOve  to  Ifae  red  dribnpiiyfi  dasoqitioa  of 

heddijr  v^eiatioa  and  repfcaemed  one  of  tiw  moat  tei^poitaftt  bands  for  v^etation  antifsis.  It  also 
provided  soil  and  geciogic  boondatf  tufomatieo.  tbt  0.49  pm  cutoff  wm  si^ificam  becaime  it 
represented  ibe  be^imtiog  of  a  apecttal  ratilge  Oom  0.£8  to  0.7^  pm  vbae  vegmatkm  reflectance 
crossovm  occur.  Generally,  this  nuge  tedpees  the  acoiCKf  of  vegetwive  investigations.  Band  4  (0.76  - 
0.90  mm  (reflecled-infiraredO  was  tuehti  for  v^etahn  hitapreiatkms  m  w^  as  soil-v^etation  contrasts. 
Band  5  (1.55  -  1.75  pm  (reflected  mid-jaflafetD  detected  tuii^dtojr  or  die  amount  of  wmer  m  v^etatkm. 
This  was  significant  for  BEXAR  because  it  indicated  {^utt  vigor  levda.  Band  5  also  provided  ^dogic 
typology.  Finally.  Band  7  (2.08  -  2.35  pm  (reflected  mid-infrared»,  which  along  with  Band  5,  had 
considerable  utility  for  the  BEXAR  analysis.  Band  7  was  the  best  TM  band  for  discriminating  geologic 
changes. 

Understanding  the  utility  of  each  band,  statistical  models  were  applied  accordingly.  The  first 
model  involved  the  testing  of  band  pte-  to  band  corrected  post-event  pair  spectral  means,  Xt^  vs. 
Xt^  respectively.  In  all,  six  band-toband  comparisons  (HO:  p|^  s  Pi^ :Hl:Pi^i«Pi^)  were  made  (i  = 

1  to  5,  and  7)  using  the  equation  below: 


where: 

a-  -  -  s  standard  deviation  of  the  variaiKe  Xl„  -  Xi  ; 
a-  and  9-.  •  sundard  deviations  of  bandi  and  band)  ; 

and 

«ad  ii|^  •  mifltiiar  of  observwions  for  banib|^  mid,^  . 

A  critical  rcgioa  based  on  equal  tail  areas  under  a  normal  curve  Xt^  -  was  chosen  for 

which  a  •  0.09.  Bfith  a  «  0.09,  then,  dm  test  sutistical  signfficanee  was  deteiained  baaed  on 
1.96*0;,^  *  “  CV,.f  ptdp,  was gie«er  titan  CV,.H0tPi^  »  Pi,,  wm  r^ecmd  and  Hl:py, 

#  Pt^  was  aoBpied  Tlie  dOfaeaeM  of  tlie  apeonl  mms  for  the  oonespaadiat  bunds,  timreliMe, 

would  be  ooiwldesed  statisifcslr  siipiiflaiu.  Hie  buMUB'^^mad  lesuf  <<  are  shown  in  TMa  . 


The  emifs  in  Tajbie  7.2.3>1  indicaitd  that  te  «a^  caae.  ^  aeaiis  of  the  band  pairinti  weie 
aigntfiniwtiy  hi  «ach  caae,  sooae  tnpe  of  ap^oal  dbante  was  detected  oatng  a  95%  ooofidenoe 

level  AJdKMih  this  testing  prooeduie  itmgrsiert  that  some  diai«e  was  appaseat,  the  amount  of  edative 
change  could  not  be  <|uantffied  using  method  Qdter  dian  piovidtag  an  iidial  owerview  of  tbe  bond 
cooqiafisons,  a  con^tariscm  of  aoeans  bjr  iisdf  is  not  coosideeed  datisticattjr  porous.  Tbe  specttal 
comparisons  did.  however,  indkate  dat  huthet  Sdiattcal  modeling  was  wvrutted 

The  nma  nmdel  used  butt  upon  die  findifl^  of  die  firai  Specffically,  dds  second  model  mnbled 
analysis  of  the  band  spectral  variance  meaaums  usif«  an  ANOVA  (analysis  of  variance).  In  the  one-way 
ANOVA  model,  the  variances  of  die  datasds  weie  laed  to  deieraaine  if  die  spectral  means  as  a  group 
were  sutisrically  differeoL  Thh  aUowed  an  assessmertt  of  the  group  diffetenoes  of  the  pre-  and  corieaed 
post-event  datasets  coduncthrely.  Such  an  aasesttnmit  was  important  because  vegeative  and  geologic 
phenomena,  as  described  earlier,  are  generally  sensed  by  several  of  the  TM  bands,  and  band  signature 
topical  overlap  needed  to  be  addressed  for  the  BEXAR  data. 

Tablo  7.2.3-1 

RMuiln  of  Tooting  tiM  Otfforoneoo  of  Bond«to-Bond  Pair  Moons 
(a  s  0.05.  n  s  6675) 


Bandsi 

cv. 

Status 

lore  vs.  IcD 

0.39 

2.01 

Reiect  HO::  Aoceot  HI: 

Zdto  vs.  2cd 

0.27 

1.21 

ReiectHO::  AocodIHI; 

3Drevs.3a» 

0.49 

2.07 

RoioctHO::  AoceotHI: 

4Drevs.4eD 

0.33 

1.39 

Roiaet  HO::  AocodI  HI : 

5wovs.5cd 

0.67 

3.74 

Retoct  HO::  Aoceot  HI: 

7pr»  vs.  7cp 

0.41 

2.40 

RsleciHO:;  AoooptHI: 

124Manl  moat  weMeompafed  oiifigd»lirpodieM»iiriiae  HOs^]^  *  14^ 

and  Hl:tti^  ^  117^ .  ibe  aest^  Sroea  die  AMOVA  tie  given  in  Table  7.23*2. 

As  indkaied  in.Tatde  7.Z5-2.  die  one-way  ANOVA  yidded  an  Kll.800^  of  11.04.  From  die 
analysts  of  variance  it  was  ii^»ied  that  the  ba«l  means,  as  a  group,  diffeied  (p<0.01).  The 
accon^ianying  statistkad  power  was  estimated  to  be  above  0.99  (00.01).  Similarly,  an  Fmax  mst  for 
honx^enehy  of  varhmoe  was  passed  at  p<0,05. 

The  ANOVA  findif^  compiimemed  die  inkial  means  testing  results.  The  F  of  11.04  was  well 
above  the  iidmied  F  of  2.25;  and  HO  was  idected,  and  HI  was  acoqited.  This  suggested  dut  dieie  were 
statistically  signifleant  spectral  differences  associated  widi  the  TM  BEXAR  dau.  The  most  compdling 
aspect  of  the  ANOVA,  however,  was  the  high  level  of  statistical  power. 

By  itself,  testing  for  group  spectral  means  did  not  provide  the  ultimate  solution  for  determining 
the  amount  and  type  of  change.  The  components  of  the  ANOVA,  though,  were  the  foundation  for  the 
final  model.  This  last  model  was  accessed  to  investigate  variant  signatures  on  a  band-by-band  matrices 
level.  This  included:  individual  band  pre-  to  band  pre-event  differences;  individual  band  pre-  to  band 
corrected  post-event  differences;  and  individual  band  corrected  post-  to  band  corrected  post-event 
differences. 

TaMg  7,2.3-2 
OiM-Wgy  ANOVA  Tabte 


'p<0.01 


The  exaoiiiuttioo  apeeanl  ^Saeeoe»  was  accomplished  usifl^  a  posteriori  IVikey  (HSD)  test 

A 

The  accompanying  yf  Qi^)  crirical  value  was  delennined  u^ig  die  ibfiowing  etpiatkMi: 


yCESD)  »  qa-oi>i.-p.,-V**®5^- 

where: 

4a-aoi.-p.«  °  ^  indexed  vahie  based  on  a  studentiaed  range  wlrere 
a  s  0.01,  p  «  die  niffiber  of  means,  and  v  -  degrees  of 
freedom  widiin  groups  Chom  ANOVA  taUe); 


MS^  =  mean  square  error  (from  ANOVA  table);  and 


n  3  number  of  observations. 


A 

With  the  BEXAR  data,  y  (HSD)  «  0.93,  p<0.01.  This  value  was  then  used  to  evaluate  the  table  of 
absolute  differences  among  means  for  the  12  bands  of  data  CTaMe  7.2.3*3). 

As  shown  in  Table  7.2.3-3,  all  but  two  diilerenoes  were  significanL  As  was  found  in  the  initial 
comparison  of  means,  the  bandi  (»e-  and  bandj  corrected  post-event  differences  were  once  again 
sutistically  significant  The  greatest  differences  among  band  pairings  were  found  in  Bands  5  0.74)  and 
Bands  7  (2.40).  Bands  1  (2.01)  and  Bands  3  C2.07)  were  also  statisticaliy  different  but  diese  were  less  than 
Bands  5  and  7.  Bands  2  (1.21)  and  Bamb  4  (1.3C0  were  significaitt,  but  these  were  only  slightly  above  the 

A 

y  (HSD)  of  0.93.  These  findings  indicaied  that  there  were  relative  vegetative  changes  from  date-to-date; 
but  the  increases  of  reflectivity,  baaed  on  Bands  1  to  4  comparisons,  indicated  that  the  changes  were  the 
result  of  normal  cyclical  growing  patterns,  .^jpaiently  die  KEXAR  events,  at  least  as  of  8  April  1991,  did 
not  damage  or  alter  the  normally-increasing  Spring  chloroplqdl  absorption  pattern  of  the  plants. 
Indeed,  the  increased  turgidity,  as  shown  by  the  Bands  5  comparison,  indicaied  that  die  vegetation  was 
entering  the  first  stage  of  a  new  annual  growth  period.  Ihis  also  was  shown  by  the  insignificant  Band  4 
pre-  vs.  Band  3  corrected  post-event  value  of  0.25  revealing  a  shift  of  healthy  vegetative  chlorophyll 
absorption  to  the  sUghtly  shorter  wavelengths. 
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Tabte7.2.3-3 

Tabto  of  AbMlute  VMurof  DNItrtnoM  Among  MMns 
MSError  «  20t.77,  n  x  6675 


Pre-Event 

B1  B2  B3  B4  B5  B7 

Corrected  Post-Event 

B1  B2  B3  B4  B5  B7 

B1 

X  40.05  26.09  23.77  9.18  27.72 

2.01  39.44  24.02  22.38  12.92  25.32 

B2 

X  14.56  16.88  «.83  12.93 

42.66  1.21  16.63  18.27  53.57  15  33 

Pre-Eveni 

B3 

X  2.32'  35.27  1.63 

28.12  13.35  2.07  3.71  39.01  0.77* 

B4 

X  32.95  3.95 

25.78  15.67  0.25*  1.39  36.69  1.55 

BS 

X  36.90 

7.17  48.62  33.20  31.56  3.74  34.50 

B7 

X 

29.73  11.72  3.70  5.34  40.64  2  40 

B1 

X  41.45  26.03  24.39  10.91  27.33 

B2 

X  15.42  17.06  52.36  14.12 

B3 

X  1.64  36.94  1.30 

Corrected 

Post-Event 

B4 

X  %.30  2.94 

B5 

X  38.24 

B7 

X 

<  V  (HSO) Of 0J3.  p<0.0l;  Th— t not ■tgnWfcart 


The  findings  in  Bands  5,  however,  indicated  another  apparent  difference.  Bands  S,  togedier  widi 
Bands  7,  suggested  tint  the  geologk:  si^puduie  of  die  area  had  been  r^eied.  By  reviewing  the  ground 
truth  informatkm  (Aiuiex  1),  this  was  detemined  to  be  tri^ered  by  a  combiiuuion  of  two  conditkms. 
First,  after  the  ea{dosion  a  large  dust  doud  was  seen  rising  up  and  moving  away  firom  die  BEXAR  site.  In 
this  cloud,  smaller  particulate  matter,  which  priev  to  the  event  was  the  surface  of  the  site,  was  removed 
After  the  event  and  the  removal  of  these  smaller  particulates,  the  surface  of  BEXAR  was  comprised  of 
larger  particulates  which  were  left  behind  Furthermore,  these  larger  particulates  had  been  moved 
around,  turned  over,  and  generally  disturbed,  changing  the  weathered-particulate  signatures.  Second 
after  the  event,  a  large  number  of  cradcs  ippezted.  on  the  surface  of  BEXAR.  These  cracks,  too,  suggested 
two  conditions:  1)  the  cracks  and  assodated  shadowing  changed  the  composite  geologic  spectral 
signature  of  BEXAR;  and  2)  the  cracking  caused  a  lateral  expansion  of  BEXAR,  and  because  of  this 
spatial  increase,  the  t|uantity  of  TM-sensiUe  geologically-related  phenomena  was  greater  in  the 
corrected  post-event  imagery  (this  expansion  concept  is  also  discussed  in  Section  7.3)- 

These  findings  suggested  that  a  combination  of  phenomenally-  and  spatially-related 
geologic'^surface  soil  changes  were  attributed  to  the  BEXAR  event;  and  other  than  natural  changes, 
minimal  (if  any)  BEXAR-related  changes  were  detected  in  the  area  vegeution.  Furthermore,  diese 
changes  were  detected  by  statistically  modeling  the  radiometric  and  spectral  resolution  characteristics 
of  the  Laruisat  V  TM  multispectral  imagery.  Itxleed,  even  with  a  relatively  small  number  of  pixel 
observations  (6675),  the  discriminating  power  of  TM  imagery  should  not  be  underestimated  from  a 
statistical  perspective.  It  was  clear,  however,  that  the  exploitation  of  these  discriminating  powers  was 
dependent  on  the  validity  and  reliability  of  the  atmospheric  removal  processes.  Thus,  within  this 
multispectral  rationale,  it  can  be  inferred  that  Landsat  V  TM  multispectral  imagery  was  successful  in 
monitoring  spectral  changes  attributaUe  to  the  BEXAR  event. 

7.2.4  AAI  Multispectral  Analyses 

There  ate  several  phenometu  associated  with  the  detonation  that  could  potentially  have 
produced  measurable  changes.  Some  of  these  were  apparent  in  the  video  recordittg  and  other  groutKl 
truth  information  provided  to  AAI  by  Autonoetric  One  was  the  creation  of  dust  plumes  at  and  around 
tlw  detonation  rite.  The  distuibanoes  which  raised  the  tfairi  may  have  (xo^ced  at  least  two  potential 
observables.  One  was  deposited  dust  in  die  vicinfty  of  and  down-wind  from  the  dust  raising  sites  (for 
example  on  vegetatioiO-  The  other  was  the  resorting  of  the  particular  siae  distributioas  at  the  detonation 
site  and  posritriy  aloi^  the  lineametu  distuibance  sites.  As  for  die  deposfted  dust,  die  resorting  probably 
enhanced  die  relative  fractions  of  fine-grained  vs.  coarse-grained  materials.  Ejqposuie  over  time  of 
mixtures  of  fine-grained  and  course-grained  materials  to  wind,  rain  and  vftxatioiid  action  generally 


wndb  to  shift  dw  more  mobUe  ftne-paiaed  fractk»s  imo  aadcs,  devices,  and  other  gnvftatkmal  wdls 
effectivdf  deoessing  dte  fraction  of  eaposed  fine-grained  matenal  per  unft  area.  Ihe  detonation  event 
toKf  have  tesiMted  die  ao&  fraedons  on  dw  qioil  pile,  enhancing  die  fraction  pm-  wdt  area  of  exposed 
fine-grained  materials. 

Such  an  enhancmnent  of  fine-grained  material  ^Ines)  on  die  spoil  pile  would  have  two 
predicted  potential  optical  effects  on  the  kugery.  First,  die  enhancemmit  of  fines  woidd  be  eiqiecied  to 
increase  the  brightness  of  the  exposed  spofi.  The  fine-pained  materials  have  higher  surfrtce  area  per 
unit  volume  dian  die  coarser-grained  mmerial,  and  the  sin^  scatterir^  (specular)  coo^ionent  erf  the 
reflected  sunlit  would  conespondirigly  be  enhanced.  Second,  die  enhanemnent  of  fines  would  diange 
the  mufttepectial  diaracteristics  of  the  material. 

The  dianges  in  muldspectral  properties  of  die  spoil  due  to  eidiancement  of  fines  are  expected 
as  a  result  of  two  factors.  One  factor  is  that  the  fines  may  be  composhionally  different,  eidier  because 
they  are  from  different  source  rodcs  or  beaujse  diey  are  more  weathered.  Different  source  materials  will 
generally  have  different  suscqxibilities  to  physical  vreadiering,  the  less  resistant  materials  produdng 
correspondingly  more  fine-grained  material  than  the  mote  resistant  source  materials.  If  the  fines  were 
derived  from  the  same  source  rocks  as  the  course  materials,  the  finer-grained  fractions  would  likely  be 
more  weathered  than  coarser-grained  fractions,  since  reduction  in  partide  siae  is  generally  a  produa  of 
weathering.  Chemical  weathering,  which  chaiiges  composition,  generally  accompanies  physical 
weathering  and  the  rate  of  diemical  weatfaeiii^  per  unit  vrrfume  is  generally  proportkxial  to  surface 
area  per  unit  volume.  Consequently  the  finer-grained  fraction  is  typically  more  chemically  weathered 
than  the  coarse-grained  fracdoa  The  rate  (rf  diemical  weathering  can  be  fast;  therefore,  even  freshly 
generated  spofl  (e.g.,  less  than  a  week  old)  can  poiemially  reveal  multispectral  effects  of  chemical 
weathering. 

The  other  factor  that  produces  diange  in  the  multispectral  properties  of  spoil  due  to  the 
enhaiKemeiu  of  the  fine-grained  fraction  is  die  change  in  effective  optical  depth-to-grain  size  ratio.  This 
changes  the  relation  of  the  single  scaoering  C^iecuiar)  and  muldpie  scattering  (diffuse)  components  of 
the  reflected  radiation.  The  diffuse  component  dominates  the  multispectral  optical  absorption 
propeities  C^iectial  structure)  of  die  maierlal,  whfle  the  specufer  conqmnent  is  dominated  by  the  index 
of  refraction  contrast  of  the  material  and  its  environment  By  enhancing  the  ffee-grain  fraction,  the 
overall  reflectance  increases  doe  to  the  increase  in  specular  lefledanoe  Qd^wr  surfitoe  area  per  unit 
vohune).  The  ^lectral  structure  wffl  inoeaae  or  deoease  dependfeig  on  wtie^er  die  optical  dqxh-to- 
gmin  size  ratio  is  closer  to  or  AMher  from  unfty  in  the  fine  grained  fracdoa 

It  is  Bkdy  that  die  spofl  mawtid  ooniaias  iron  oaide  as  an  aooeasory  pigmenting  adKfal;  since 
itisc»ieof  diemoittibitpiliBMiofwsthstlmpiotlactttandltoecaialnaMStaoila.  feonosMBlMtta  vwy 
small  optical  depth  Ctdongif  shsotbing)  m  dm  vMMe  petdon  of  dm  apstawaa^  and  dm  fine-grain 
fraction  shmdd  diamfeie  hem  enhanced  hen  oaMh  specand  Mmctuie  tcitdve  to  dm  mote  coarse- 
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gnined  CnctkMi.  If  true,  a  potential  obaervab'e  resulting  from  a  deposit  of  dust  or  a  resorting  of 
materials  following  detonatimi  would  be  a  spectral  'reddening*  of  the  spoil  material.  The  spoil  after 
detonation  would  not  rmly  be  bri^ter.  but  tbe  irtfnued  bands  would  increase  in  reflectance  more  than 
the  visible  bands  in  a  manner  consistent  with  the  spectral  structure  of  iron  oxide.  The  presence  of  other 
materials  with  spectral  structure  may  complicate  the  chartge;  however,  as  for  most  soils  that  contain  low 
organic  content,  the  spectral  "reddening*  of  the  spoU  should  be  the  principal  observable  resulting  from 
an  enhancement  of  the  fine-grained  fraction,  if  there  was  one. 

In  addition  to  changes  in  the  spoil  pile,  the  area  down-witxl  of  the  spoil  pile  and  in  the  vicinity 
of  the  dust-raising  lineaments  may  also  have  changed.  The  principal  observaUe  would  be  the 
modification  of  the  background  spectral  properties  by  the  addition  of  d^x»ited  dust  Although  such  a 
change  is  potentially  observable,  the  natural  changes  in  the  vegetation  spectral  properties,  associated 
with  changes  in  vegetation  density  and  growth  cycle,  and  in  the  relative  fractions  of  illumiruted 
vegetation  and  unvegeuted  floor  rruterial  due  to  changes  in  sun  angle,  nuke  detection  of  change  due  to 
added  dust  unlikely.  Spectral  'reddening*  of  the  spoil  pile  is  a  substantially  less  ambiguous  potential 
observable  related  to  the  prcxluction  of  dust. 

Another  (henomenon  that  may  have  prtxluced  an  observable  was  vegetal  stress  induced  by  the 
disruption  of  the  rcxx  systems  during  the  event  There  are  three  factors  that  militate  against  use  of  that 
potential  observable  in  this  study,  however.  First,  most  of  the  area  directly  affected  by  the  detonation 
was  unvegeuted  spoil  material.  Although  there  may  have  been  significant  disruption  of  the  surroutKiing 
vegeution,  the  area  of  greatest  disruption  was  the  spoil  pile.  Whether  there  was  significant  disruption  of 
the  root  systems  at  distances  beyond  the  spoil  pile  is  uncertain.  Second,  changes  induced  by  stress  can 
not  unambiguously  be  isolated  from  changes  induced  by  natural  effects  unrelated  to  the  event  Third, 
the  four-day  time  pericxl  from  the  evem  to  the  collection  of  the  post-event  TM  image  was  shon  in 
comparison  to  the  probably  longer  time  scales  of  stress-induced  spectral  change  (exclusive  of  severing 
the  water  supply). 

A  third  potential  (4>servable  is  measurable  expansion  of  the  spoil  pile  produced  by  the 
detorution.  Edge  expansion  could  have  potentially  CKCurred  as  a  result  of  subsidence  or  ballistic 
expansion.  The  subpbcel  detection  capability  of  AASAP  potentially  makes  possible  the  detection  of  such 
change  at  spatial  scales  that  are  smaller  than  conventional  whole  pixel  techniques,  atxi  this  potential 
observable  is  discussed  in  mote  detail  in  Section  7.3 

Other  potential  observables  may  have  resulted  from  changes  in  fracturing  of  the  spoil  pile  and 
a  variety  of  other  possible  effects  discussed  above.  The  two  principal  observables  discussed  here, 
enhancement  of  spoil  pile  fines  and  spoil  pile  edge  expansion,  were  the  focus  of  the  subcontraaor's 
search  for  diange.  Possible  evidenoe  of  other  chenge  was  conskieted  secondarily. 

To  search  for  possible  evidenoe  of  change  induced  by  dqxwhed  dust  or  mcierial  tesbrting.  the 
propoMd  spoU  "reddening*  obsmvable  was  souf^  hi  the  data.  The  ground  tnith  photographic  dau  were 
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used  to  identiiy  pixds  estimated  to  be  at  or  near  dte  detorouion  ground  zero  site.  'Ibiee  pizeb  were 
sdected  in  the  pie-eveitt  image,  and  the  identical  (widiin  registration  accuracy)  pixds  were  selected  in 
the  po«-event  image.  The  two  sets  of  piasels  are  lismd  in  TaWe  7.2.4-1.  The  three  enviroimierttal 
correoions  KMed  in  Tadiles  7.2.2-1,  7.21-2,  and  7.2.2-4  were  applied  to  die  pixels  in  TaWe  7.2.4-1.  Each 
corrected  piiml  in  the  post-event  scene  was  divided  by  the  corresponding  pixel  in  the  pre-event  scene  to 
derive  a  pixel  ratio  (post/pte)  for  each  of  the  duee  selected  (^xel  locations.  These  are  shown  in  Table 
7.2.4-2.  Also  included  is  the  average  of  the  three  i^xd  ratios.  With  the  corrections  that  were  applied, 
these  pixel  ratios  are  proportional  to  the  ratios  of  the  image  i^ane  reflectances  of  the  after  material 
relative  to  the  before  material.  The  results  are  generally  consistent  with  enhanced  reactance  in  the 
infrared  bands  in  the  post-  relative  to  the  pre-  pixels,  compared  to  the  ratios  in  the  visitde  bands  where 
reflectances  are  mote  comparable. 

Tai»l«  7.2.4-1 
FIhm  EntiancMmiM  Pht«l« 

(Data  Humbars) 


Landsat  Thenutic  Mapper 

Imann  Plana 

7  Match  1991 

Ptw»la/nM\ 

8  April  1991 
HMlS.tDN) 

1 

117 

126 

135 

150 

149 

145 

2 

62 

67 

72 

84 

82 

80 

3 

97 

106 

110 

132 

129 

127 

4 

81 

94 

99 

111 

113 

113 

5 

119 

149 

149 

186 

191 

192 

7 

77 

89 

88 

111 

118 

117 

Tabla  7.2.4-2 

Pixal  Rmoa  Aflar  Conaetiona 
(Plnaa  Enhancamant) 


Landsal  ThamaHc  Mapper 

Imarw  Plano 

Aftar/Befora 

Bi»aian>wt 

Averaoa 

PhtaiRmiQa 

1 

1.13 

.99 

.85 

.99 

2 

1.16 

1.02 

.90 

1.03 

3 

1.14 

1.01 

.95 

1.03 

4 

1.15 

1.00 

.94 

1.03 

5 

1.29 

1.05 

1.08 

1.14 

7 

1.18 

1.09 

1.09 

1.12 

» 
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The  results  are  further  consistent  with  imagery  of  the  small  poods  near  the  detonation  site.  The 
(tetonation  wtHild  be  expected  to  disturb  the  floor  material  in  the  ponds,  aitd  as  the  material  settled  the 
fines  would  settle  out  later  than  the  coarse-ground  material.  Some  of  the  silt  likely  remained  suspended 
in  the  water  fdlowing  the  detonation.  The  result  should  be  a  fines  enhancement,  and  the  pixel  ratio 
(post/pre)  should  similarly  show  the  proposed  charar^ristic  ’reddening”  observaUe.  Pbcels  from  the 
ponds  were  selected  from  the  pre-  and  post-event  images,  arxl  they  are  listed  in  Table  7.2.4-3.  The  pixels 
were  environmeiually  corrected  and  post-/pre-event  pixel  ratios  were  calculated  directly  analogous  to 
the  processing  of  the  spoil  pixels  in  Table  7.2.4-1  that  led  to  the  results  in  Table  7.2A-2.  The 
corresponding  results  for  the  spoil  pond  pixels  are  shown  in  Table  7.2. 4-4.  The  reddening  was  stronger, 
consistent  with  the  reduced  index  of  refraction  contrast  within  the  water.  The  relatively  high  Band  7 
value  also  suggested  that  the  silt  may  have  been  partially  suspended. 


TabI*  7.2.4-3 

Spoil  Pond  Fines  Entiancament  Pixels 
(Data  Number) 


Landsat  Thematic  Mapper 
Imaoe  Plane 


7  March  1991 
Pixels  fPNl 


8  April  1991 
Pixels  fPNl 


1 

70 

79 

88 

97 

96 

92 

2 

47 

52 

57 

67 

65 

63 

3 

86 

95 

99 

120 

117 

115 

4 

76 

89 

94 

106 

108 

108 

5 

118 

148 

148 

185 

190 

195 

7 

76 

88 

87 

110 

117 

116 

Table  7.2.4-4 
Pbcal  Ratios 

(Pond  Fines  Enhancement) 


Landsat  Thematic  Mapper 

After/Before 

Average 

Imaoe  Plane 

Pixel  Ratios 

Pixel  Ratios 

1 

1.01 

1.10 

1.19 

1.10 

2 

1.10 

1.12 

1.15 

1.12 

3 

1.11 

1.21 

1.15 

1.12 

4 

1.24 

1.04 

1.15 

1.15 

5 

1.16 

1.78 

1.53 

1.49  ^ 

7 

1.50 

1.85 

1.55 

1.63 

i 
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The  findings  appeared  to  support  the  {Mt^josal  that  fines  enhancement  nay  have  been  the 
source  of  the  otoerved  reddening.  In  particular,  die  finding  that  the  bright  spot  (indicative  of  a  local 
concemration  of  a  fine-grained  materiaO  was  redder  than  the  avenge  spoil  ^{^>oned  the  {m^xxal  that 
the  fines  were  reddm’  than  the  average  spcul.  Ibe  observation  that  the  reddening  was  enhanced  in  the 
pond  relative  to  da  open  spoil  was  consistent  with  the  expected  enhancement  of  the  difiiise  component 
relative  to  the  specular  component  of  scattered  mdiation  in  water  vs.  air,  since  the  water  suppresses  the 
index  of  refraction  contrast  between  the  partides  and  their  environment  and  hence  suppresses  da 
specular  component  of  scattered  radiation. 

Additional  effort  is  needed  to  verify  whether  the  (diservable  is  a  material  {dienomenon 
associated  with  the  detonation  event,  but  the  results  are  at  least  consisant  with  the  proposal  that  such  an 
observable  may  be  a  potential  indicator  phenomenoa 
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73  AAI  MUUnSPECTRAL  SUBPDCEL  ANALYSES 


As  mentioned  previously  in  SectkMis  5.3  and  7.2.4,  the  BEXAR  detonation  may  have  caused 
changes  in  the  nature  and  distribution  of  surface  material  aloi^  the  edge  the  site’s  qxxl  {^e  that 
could  be  detected  using  subpixel  technkpies.  For  instance,  a  pixel  at  the  edge  of  the  spc^  pile  mi^t  be 
located  so  that  the  spoil  covers  half  of  the  pixel  area.  The  remaining  pixel  area  mi^  be  coveted  by 
many  non-spoil  surface  materials  includit^  vegetatkm,  dirt,  and  rocks.  A  pixel  sudi  as  this  was  defined 
here  as  a  mixed  pixel.  It  contained  both  subpixel  quantities  of  a  material  interest  or  "target*  On  this 
case  spoil)  and  other  urget  background  nuterials.  Quantifying  the  subpixel  and  the  whcde  pixel  target 
amounts  using  the  pre-  to  post-event  TM  data,  an  examination  of  the  spoU  pile  edge  expansion  was 
conducted. 

To  search  for  this  subpixel  change,  the  AASAP  was  tised.  AASAP  provides  information  about  the 
composition  of  mixed  pixels  that  can  be  used  to  identify  and  classify  change.  AAI  originally  designed 
AASAP  as  an  automatic  target  detection  algorithm  for  multispectral  image  processing.  AASAP  prcxesses 
an  image  on  a  pixel-by-pixel  basis  to  identify  the  location  of  specific  materials  or  targets  of  interest  For 
a  given  pixel,  AASAP  first  estimates  the  spectral  contribution  of  non-urget  materials  (background)  and 
then  removes  the  backgrouixl  contribution  to  reveal  the  presence  of  target-like  materials. 

One  output  of  the  AASAP  background  removal  process  is  the  K-faaor.  The  K-factor  represents 
the  fraction  of  badcgrouixl  removed.  It  provides  a  measure  of  the  area  of  a  pixel  covered  by  background 
material  but  is  more  properly  defined  as  the  fraction  of  radiance  contributed  by  background  in  the 
original  pixel  spectrum;  the  remaining  fraction  (1-k)  is  radiance  contribution  of  the  target  material. 

A  subpixel  change  detection  technique  using  the  Landsat  TM  dau  is  described  in  this  section. 
The  AASAP  K-faaor  output  was  used  to  show  that  for  selea  jmcls  located  on  the  spoil  pile  edge,  the 
subpixel  amounts  of  spoil  material  increased  in  the  after  image  relative  to  the  before  image.  The 
subpixel  occurrence  of  more  spoil  material  in  the  after  spoil  pile  edge  pixels  indicated  the  spoil  pile 
expanded  after  the  detonation. 

The  first  step  in  this  technique  was  to  define  which  pixels  were  located  on  the  spoil  pile  edge. 
Figure  7.3-1  is  a  Landsat  TM  image  with  two  spoil  piles  identified.  Spoil  pile  *1  (BEXAR)  is  near  the 
center  of  the  image.  Spoil  pile  #2  is  in  the  southwest  comer.  Spoil  pile  *2  was  from  a  previous  test  site, 
and  it  was  assumed  no  dimensional  change  had  occurred  from  the  pre-  to  post-event  TM  collections. 
On  and  around  each  of  the  two  spoil  piles,  several  16-pixel  windows  were  selected  as  sam{^  to  develop 
the  spectral  dau  and  test  for  spoil-pile  expansion  (Figure  7.3-2). 
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Figure  7.3-2 

Pixel  Samples  On/Around  Spoil  File  Sites 


These  patdi  positions  were  setecied  using  visual  determination  of  the  ty^jraximate  spoil  edge 
location  from  disphiys  of  the  six  TM  images.  Six  of  the  windows  were  placed  over  borizontally-oiiented 
edges  in  the  image,  and  three  patches  were  pfaioed  over  veiticaBy-oriented  edges. 

The  nnt  step  was  to  idmttify  sin^  pixel  locations  along  the  qwil  pile  edges  within  eadt  row  (tf 
the  windows  (or  each  column  in  the  case  of  die  three  windows  on  vertical  edges).  Each  tow  (or  column) 
in  the  windows  represented  a  sin^  pixel  transect  aaoss  the  edgt.  This  was  accomplished  using  AASAP 
K-£actor  analyses.  Again,  the  target*  of  interest  in  this  case  was  the  spoil  material.  The  K-Eact(» 
represents  die  fraction  of  non-spod  material  m  die  pixd  being  tested.  A  large  K-Eicror,  for  instance  0.70, 
means  a  large  amouiu  of  background  Oipproxurnttely  7096  of  the  spectral  energy  from  the  f^xeD  is 
present  in  the  pixel.  K-^ctors  are  provided  in  Figure  7.3-3  for  die  windows  on  spoil  pile  #1  from  the  pre- 
event  image  (if  the  background  material  could  not  be  identified  a  dash  was  entered  in  the 
corresponding  pixeO.  Note  the  gradual  change  in  the  K-fictor  from  the  inner  pixels  to  the  carter  pixels. 
For  additional  disoBsion,  the  *yellow*  window  (top  lefD  frcmi  spoil  pile  #  1  is  rf^raifavi  in  Figure  7.3-4. 

In  the  second,  third,  and  fourth  row  of  the  patch,  the  K-factors  in  die  right  three  columns 
gradually  increased  from  the  inside  ctf  the  spoil  pik:  to  the  outside.  Moving  to  the  left,  the  K-friaors  first 
increased,  indicating  that  mote  and  more  background  material  had  to  be  removed  to  indicate  the 
presence  of  the  spoil  material  This  gradual  change  in  die  K-factor  was  one  indication  that  the  spoil 
material  (or  spoil-like  materiaO  transitioned  gradually  into  the  background  At  this  location  it  was  not  a 
sharp,  distinctive  spoil/badcgtound  edge.  In  rows  3  and  4,  the  K-friaor  peaked  and  then  decreased  This 
indicated  the  location  of  an  edge.  The  maximum  K-factor  alcxig  each  transect  is  identified  in  Figure  7.3- 
5  by  cross  hatches. 
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Ftgur»  7^.  K  Factor  Analysis  (Betopo  Detonation) 
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FIgurt  7.341.  K  Factor  Analysis  (Before  Dtlonation) 


FlfHI*  7.3-5 

k-fmm 

r  Amlytn  (F—t  Dtowtlon) 

The  0.29  K-fiictor  in  the  boctom  teft  was  associated  mdi  a  relatively  poor  detection  confidence 
level.  The  confidence  level  is  determined  as  an  output  parameter  (rf  AASAP.  This  low  detection 
confidence  value  corresponded  wtth  a  low  omifidence  that  spoQ  nuUerial  was  present  The  low  detection 
ccMifidenoe  may  have  indicated  the  presence  of  dirt  or  some  other  material  that  was  target-lifee,  but  not 
spoil,  and  therefore  could  be  overiooked. 

Spoil  pile  expansion  may  produce  a  "subpixel*  change,  where  subpixel  change  is  defined  as  a 
change  in  the  mixture  of  materials  that  combiiw  to  produce  the  multispectral  dau  nundrers  for  a 
specific  pixel  Thus,  if  the  edge  of  the  spoil  pile  expands,  as  shown  in  Figure  7.3^,  then  the  fraction  of 
spoil  in  the  pixel  will  increase  as  illustrated  in  Figures  7.3-6b  and  7.3-6c  The  AASAP  K-factors  (xovide  a 
direa  measure  of  the  amount  of  a  material  of  imerest  in  the  pixel  under  evaluation. 

The  K-factw  edge  detection  tedmique  was  used  on  the  top  three  pixel  wukIows  shown  in  Figure 
7.3-2  to  identify  spoil  pile  edge  expansion  near  die  drill  hole  in  spoil  pile  *1.  As  a  control,  the  three 
pixel  windows  at  the  bottom  of  spoil  pile  *1  and  the  three  pixel  wirxiows  at  spoil  pile  *2  were  used  to 
demonstrate  that  the  spoil  edge  did  not  expand  in  these  control  locations. 

Figure  7.3-7  provides  the  K-£actor  values  for  the  six  pixel  windows  for  spoil  pile  *1  after  the 
detonation,  with  the  maximum  K-factors  indicated  by  cross  hatches.  Figures  7.3-fi  and  7.3-9  provide  the 
K-factor  values  for  the  pixel  windows  from  spoil  s2  before  and  detonation,  respectively. 


Flguft7^.  Spoil  Expansion  (Detail) 


Plgyra  7.»« 

SpoN  ijqMiitlen  (OMM) 


These  figures  represented  the  results  of  the  spoil  pile  edge  expansion  analysis.  By  comparing 
the  position  of  tte  spoil  edge  in  the  pre-event  image  (the  hatched  pixels  in  Figure  7.3-5)  to  its  position 
in  the  post-event  image  (the  hatched  pixeb  in  Figure  7.3-7),  it  was  concluded  that  the  spoil  pile  edge 
expanded  outwardly.  The  pixel  locations  of  the  maximum  K-faaors  in  the  top  right  and  top  left  pixel 
windows  for  spoil  pile  #1  were  shifted  outwardly  in  the  post-event  image.  This  indicated  an  erqiansion  of 
the  spoil  pile  near  the  detonation  site.  No  consistent  evidence  of  expansion  was  apparent  in  the 
remaining  control  pixel  patches  at  the  bottom  of  spoil  pile  #1  and  at  spoil  pile  *2  (Figures  7.3-5,  7.3-7, 
7.3-8,  and  7.3-9).  This  supported  the  hypothesis  that  the  spoil  pile  expanded  near  the  detonation  site. 
Figure  7.3-10  provides  a  schematic  of  spoil  pile  *1  with  the  pixel  shown  in  the  Figure  7.3-6  (a  to  c) 
indicated.  The  expansion  near  the  top  of  the  spoil  is  shown  schematically  and  does  not  correspond  to 
the  actual  degree  of  exfumsion. 


) 
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Figure  7.3-10.  Spoil  Expansion 


7.4  HYPERSPECTRAL  ANALYSIS 


Th  hyperspectral  component  of  the  lesearch  was  accomplished  with  data  collected  using  the 
MTL  Systems,  Inc  63-band  ASIS.  As  a  subcontractor  to  Autometric,  MTL  flew  the  ASIS  several  times 
during  the  first  week  of  April  1991.  Specifically,  three  fli^us  were  made  befme  the  event  on  2  April,  and 
one  was  made  after  the  event  on  4  April. 

The  63-band  ASIS,  flown  in  a  twin-engined  Piper  Aztec  airplane,  collects  l6-bit  etectromagnetic 
energy  in  the  0.4  to  2.5  |im  range.  Three  grating  spectrometers  operate  in  the  ASIS.  One  (a  silicon 
detector)  provides  24  bands  of  data  in  the  0.4  to  1.0  pm  spectral  region,  and  the  band  wi^  for  these  24 
is  approximately  25  tun.  The  second  spectrometer  covers  the  1.1  to  1.8  pm  region  with  seven  bands,  and 
each  of  these  has  an  approximate  spectral  width  of  125  nm.  The  bandwidth  of  this  spectrometer  is  wider 
to  compensate  for  lower  available  sigtul  inputs  resulting  from  atmospheric  absorptions.  The  third 
spectrometer  provides  32  bands  in  the  2.0  to  2.5  pm  tegioa  Each  of  these  has  an  approximate  spectral 
width  of  16  tun.  The  second  and  third  spectrometers  use  lead  sulfide  detectors  and  are  cooled  by  two- 
stage  thermoelectric  cooling  devices. 

Summaries  of  the  four  ASIS  collections  of  BEXAR  are  {Movided  in  Table  7.4-1. 


Table  7.4-1 

Airborne  SpcctrtMndlonietric  Imaging  System  (ASIS) 
BEXAE  Image  CoUeetton  Summaries 


Date 

Time 

MTL-derived 

Altitude 

Spatial 

2  April  1991 

11:30 

BEXAR  2 

13,000 

6  by  6  m* 

2  April  1991 

11:40 

BEXAR  4 

13,000 

6by6  m 

2  April  1991 

11:55 

BEXAR  1 

17,000 

10  by  10  m 

4  April  1991 

14:45 

BEXAR  5 

13,000 

6  by  6  m* 

*Sent  to  AAI  (see  Section  7.4.2) 

As  with  the  Landsat  TM  imagery,  both  Autometric  and  AAI  analyzed  the  ASIS  data.  Accordingly, 
these  analyses  have  been  repotted  separately  and  are  provided  in  the  following  Sections  7.4.1  and  7.4.2, 
respectively. 


7.4.1  Autometric  Hyperspectral  Analyses 


Airborne  spectroscopic  devices  were  originally  designed  to  satisfy  discriminate  analytical 
requirements.  The  narrow-band  discrimination  capability  provided  by  the  systems  compliment 
statistic  ’.ly-driven  multispectral  imagery  (e.g.,  Multispectral  Scanner  and  TM)  methodologies. 
Discriminate  techniques  (mainly  spectral  response  curve  investigations)  are  commonly  used  by  analysts 
reviewing  the  voluminous  data;  indeed,  statistical  models  built  with  50  to  200+  spectroscopic  bands 
(input  variables)  are  difficult  to  apply.  Accordingly,  for  this  portion  of  the  research,  discriminate 
spectral  curve  analyses  were  performed  on  pre-  and  post-event  ASIS  hyperspectral  data  for  studying 
event-related  changes. 

Before  the  spectral  daot  could  be  compared,  several  preparatory  steps  were  required.  For 
instance,  the  originai  ASIS  imagety  was  collected  as  l6-bit  data.  The  display  devices  at  the  Autometric 
image  processiitg  laboratory,  however,  being  standard  8-bit  black  and  white  and  24-  and  32-bit  color 
displays,  were  not  compatible  with  the  raw  l6-bit  data.  The  first  step  in  the  ASIS  processing,  therefore, 
was  a  reduction  of  the  l6-bit  (0  to  65535)  imagery  to  8-bit  (0-255).  This  was  performed  separately  on 
each  band  using  the  band-uniqud  statistical  profile  information.  For  a  given  l6-bit  band,  minimum  ?nd 
maximum  values  were  declared  at  2o  about  the  mean  (i.e.,  approximately  95%  of  the  data).  These  were 
divided  by  256,  and  the  resulting  values  became  the  absolute  minimum  and  maximum  for  the  output 
band.  The  remaining  pixels,  starting  with  the  declared  minimum  value,  were  density  sliced  into 
brightness  level  interval  classes,  each  having  a  l6-bit  range  of  256.  The  last  class  was  comprised  of  a 
remainder  amount.  In  turn,  these  interval  classes  were  used  to  remap  each  l6-bit  pixel  to  an  8-bit  scale. 
Although  this  process  reduced  the  original  radiometric  range,  the  output  8-bit  bands  nevertheless 
retained  the  overall  l6-bit  histogram  profile  as  well  as  vital  relative  band  spectral  integrity.  This 
conversion  was  performed  for  each  of  the  four  63-barxi  ASIS  datasets. 

After  being  converted  to  8-bit  imagery,  the  four  ASIS  sets  were  reviewed  for  image  qualities.  An 
example  of  Band  10  for  each  collection  is  shown  in  Figure  7.4.  l-I.  (BEXAR  3  was  a  null  scene  due  to 
sensor  malfunctions,  hence  the  numbering  BEXAR  1,  2,  4,  and  5.)  Specifically,  keeping  in  mind  the  goal 
of  this  analysis  was  to  detea  spectral  changes  caused  by  the  BEXAR  event,  these  four  datasets  were 
analyzed  to  determine  which  two  would  be  most  appropriate  for  making  pre-to-post-event 
discriminations.  There  was  no  choice  for  the  after  image;  BEXAR  5  was  the  only  post-event  image.  For 
the  pre-event  scene,  though,  there  were  three  potential  choices.  BEXAR  4  was  eliminated  because  a  large 
cloud  shadow  obscured  the  general  test  site  area.  Considering  BEXAR  5  had  a  spatial  resolution  of  6m,  it 
was  decided  that  BEXAR  2  —  also  at  6m  and  approximately  the  same  geometric  orientation  as  BEXAR  5 
—  would  serve  as  the  pre-event  dau.  BEXAR  1,  at  10m  and  skewed  badly,  would  have  required 
substantial  rectification,  and  this  would  have  signiftcantiy  degraded  the  spectral  integrity  of  the  dau 


77 


As  shown  in  Figure  7.4.1-1,  BEXAR  2  did  contain  a  series  of  scanline  dropouts,  but  these  were  not  over 
the  BEXAR  site.  Consequendy,  they  were  ignored. 

Once  the  pre-  (BEXAR  2)  and  post-event  (BEXAR  5)  dau  were  chosen,  the  image  sets  were 
registered.  This  ensured  for  a  given  ground  location  (X,Y)  in  the  pre-event  image,  the  same  specific 
ground  location  (X,Y)  could  be  located  in  the  post-event  inuge.  The  specific  method  of  registration  was 
discussed  in  Section  7.1;  with  the  ASIS  data,  however,  the  pre-event-  imagery  served  as  the  geometric 
ground  truth  (instead  of  a  map  as  in  the  multisensor  registrations),  bi  turn,  the  post-event  imagery  was 
registered  to  the  pre-event  base  image.  To  achieve  this,  12  GC3*s  were  selected  from  the  general  test  area, 
and  using  them  in  the  rectification  model,  an  RMS  error  of  <1.50  pixels  was  obtained. 

In  reviewing  the  65  bands  of  each  dau  set,  it  became  apparent  that  the  ASIS  sensor  gains  were 
set  differently  from  coilection-to-coUection.  Upon  follow-up,  the  MTL  sensor  operator  reported  these 
different  gain  settings  were  an  effort  to  boost  the  input  signal  during  the  post-event  collection. 
Apparently,,  this  was  done  to  compensate  for  thin  upper-atmo^heric  clouds  on  the  afternoon  of  4  April 
which  absorbed  significant  portions  of  near-  and  mid-infrared  energy.  Thus,  the  post-event  dau  were 
much  bri^ter  relative  to  the  pre-event.  To  remove  these  impacts  from  sensor-related  differences,  two 
calibration  sites  were  chosen  away  from  the  immediate  BEXAR  site.  One  was  a  vegeution  sample,  and 
the  second  was  a  bare  soil  sample.  Gain  offsets  were  then  determined  for  each  band,  and  these  values 
were  used  to  calibrate  the  post-event  dau  to  a  pre-event  baseline. 

Results  from  the  data  rectification  and  calibration  preparatory  functions  are  shown  in  Figure 

7.4. 1- 2  using  another  Band  10  example. 

To  analyze  the  discriminate  sigtuture  features  of  the  pre-  and  post-event  ASIS  daU,  nine  ground 
locations  were  determined  from  the  registered  imagery.  Samples  consisting  of  both  ’soil”  and 
"vegetative"  areas  were  uken  from  the  BEXAR  site  and  surrounding  areas.  These  are  shown  in  Figure 

7.4. 1- 3.  One  of  the  test  samples  (site  6)  was  under  the  leading  edge  of  a  cloud  shadow  on  the  pre-event 
imagery  and  was  not  included.  Hence,  a  toul  of  eight  valid  samples  as  used  for  the  discriminate 
sigruture  evaluations. 


Figure  7.4,1 -2 
Rectlfisc!  and  Calibrated 


Airborne  Spectroradlometric  Imaging  System  (ASIS)  Band  10  Examples 
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Figure  7.4.1-3 

Soli  and  Vegetation  Sample  Areas 
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Complementing  previous  study  results,  an  analysis  of  the  soil  and  vegetative  sample  spectral 
curves  from  the  two  collections  revealed  the  detection  of  apparent  event-related  and  event-unrelated 
changes.  Unfortunately,  the  spectra  suggested  the  ASIS  sensor  was  not  functioning  at  full  capacity. 
Specifically,  in  ail  sam{de  cases  Bands  32  to  63  had  extremely  low  brightness  value  counts;  they  were 
also  highly  correlated  (>95%).  Bands  32  to  63  were  collected  by  the  third  spectrometer,  and  rHeir 
corresponding  minimal  utility  indicated:  1)  spectrometer  three  was  not  functioning  properly;  2) 
spectrometer  three  was  functioning  as  designed,  but  it  was  not  sensitive  enough  to  detect  the  variant 
spectral  information  resident  in  the  2.0  to  2.S  pm  spectral  region;  or  3)  the  upper-atmosphere  clouds 
absorbed  too  much  of  the  near-  and  mid-infrared  energy,  and  not  enough  energy  was  reaching  the 
sensor.  Consequently,  considerable  amounts  of  data  on  the  geologic  surface  composition  were  not 
available.  Comparatively,  ASIS  spectrometers  oocuind  two  apparently  functioned  normally.  They 
produced  valuable  dau  from  0.4  to  1.0  pm  and,  mote  importantly,  in  the  1.1  to  1.8  pm  range. 

With  regard  to  the  detection  of  diange,  two  example  plots  are  shown  in  Figures  7.4. 1-4  (site  7  - 
vegeution)  and  lA.\-‘i  (site  5  -  soiO-  Figure  7.4.1-4,  together  with  other  plots  not  shown,  indicated  no 
apparent  change  could  be  detected  from  0.4  to  1.0  pm.  This  suggested  the  chlorophyll  absorption  levels 
remained  constant  from  the  2  to  4  April  collection  dates.  Slight  spectral  changes  were  observed,  but 
these  were  not  considered  significant.  In  contrast,  change  in  the  1.1  to  1.8  pm  portion  was  noticed. 
Specifically,  a  slight  increase  in  the  corrected  pc»t-event  dau  indicated  a  small  increase  in  turgidity. 
The  greater  turgidity  level  was  likely  the  combined  result  of  a  continuing  seasonal  natural  increase  and 
also  the  later  time-of-day  sensing  on  4  April.  Because  of  the  lack  of  dau  in  the  2.0  to  2.5  pm  range, 
however,  no  reliable  assessment  of  the  spectral  impact  of  the  detonation  shock  dust  cloud  could  be 
made.  Indeed,  some  inaease  in  reflectivity  frran  the  post-event  collection  in  the  1.1  to  1.8  pm  range 
could  have  been  attributed  to  particulate  deposits  on  the  vegeution  (especially  for  sites  7  and  8),  but 
none  was  detected.  (Given  the  high  wind  velocity  at  the  site,  the  dust  coating  had  probably  been  blown 
off  the  surface  of  the  evergreen  vegeution  soon  after  the  test.) 

Although  natural  change  was  noticed,  no  event-related  changes  were  detected  using  the 
vegeution  sample  sites.  These  results  agree  with  the  findings  from  the  multispecual  analyses  (Section 
7.2). 

For  the  soil  sites,  the  dau  seemed  to  contradict  results  from  previous  analyses.  In  particular, 
post-event  brightness  values  in  the  1.1  to  1.8  pm  range  were  less  than  values  in  the  pre-event  imagery. 
Again,  though,  the  2.0  to  2.5  um  range  was  not  usable  for  comparison.  Referring  to  Sections  7.2  and  7.3 
and  Annex  1,  however,  three  physical  changes  to  the  test  site  were  detected  previously:  1)  expansion  of 
the  test  site  surface  (as  shown  in  TM  Band  7  increases  in  the  ANOVA  and  Tukey  HSD  tests  atKl  the 
subpixel  processing);  2)  the  alteration  of  the  surface  particulate  material  (as  shown  by  the  spectral 
reddening  in  TM  Bands  5  and  7);  and  3)  the  formation  of  a  .series  of  .surface  cracks  after  the  explosion. 
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Figure  7.4.1«4.  Site  7  -  Vegetation 


The  ASIS  data  had  a  spatial  resolution  of  6m,  and  this  was  mudi  hi^r  than  the  30m  of  the  TM 
imagery.  Compared  to  the  TM  imagery,  with  this  increased  resolution  the  cracks  and  accompanying 
shadow  features  comprised  more  of  an  ASIS  composite  pixd  signature.  In  the  6m  ASIS  imagery,  these 
changed  sutes  reduced  considerably  the  level  of  sensed  energy  during  the  post-evetu  collection. 
Furthermore,  the  shadowing  effea  was  also  amplified  by  the  low  sun  position  during  the  14:45  Qocai) 
collection  time. 

Based  on  these  results,  it  was  inferred  that  data  from  the  ASIS  did  detect  specific  changes  related 
to  the  underground  detonation.  The  greater  spectral  resolution  of  the  ASIS  data,  however,  did  not 
provide  substantially  more  information  compared  to  the  six-barui  Larxlsat  TM  imagery.  This  may  not 
have  been  the  case  had  spectrometer  three  provided  usable  data.  Also,  some  information  was  lost 
during  the  l6-bit  to  B-bit  conversions,  as  well  as  the  rectifications.  In  contrast,  the  considetabie  spatial 
resolution  increase  did  detect  the  cracking  and  assoaamd  shadowing,  and  these  features  were  not 
isc^ated  using  the  TM  data. 

Perhaps  additional  processing  methodologies  should  be  considered  to  utilize  fully  an  ASIS  or 
other  spectroscopic  data  set.  Indeed,  a  combination  of  16-^1  stochastic  time-  and  frequency-domain 
modeling  techniques  with  conjunctive  8-bit  discriminate  arid  limited  statistical  image  processings  may 
provide  a  more  rigorous  analyticai  approach.  In  general,  though,  the  results  from  the  8-bit  discriminate 
methods  used  here  agreed  with  the  findings  from  the  other  components  of  the  research:  spectral 
changes  related  to  the  underground  detonation  were  detected  using  remotely-sensed  dau  together  with 
ground  truth  data. 


1A2  AAI  Hypeispectral  Analyses 


Complimenting  the  multispectral  (Section  7.2)  and  subpixel  (Section  7.3)  analyses,  AAI  sought 
another  line  of  evidence  for  potential  fines  enhancement  using  the  ASIS  data.  As  indicated,  the  ASIS 
dau  had  higher  spectral  and  spatial  resolutions  than  the  Landsat  imagery. 

Not  all  of  the  ASIS  data  was  found  to  be  usable  as  illustrated  in  Figure  7. 4.2-1.  In  the  figure  the 
data  number  ranges  for  the  pre-event  2  April  1991  and  post-event  4  April  1991  images  are  shown. 
Unfortunately,  many  of  the  bands  from  spectrometer  could  not  be  used 

In  Figure  7. 4.2-2  are  shown  the  averages  of  18  pixels  (6x3  grids)  spectra  for  a  spot  on  an 
unvegetated  road  area  away  and  upwind  from  the  spoil  pile.  Two  sets  of  18  pixels  were  obtained,  one  set 
from  the  pre-event  image  and  the  other  from  the  post-event  image.  It  was  attempted  to  extract  the  pixels 
from  the  identical  locations  on  the  road  in  the  two  images.  It  was  apparent  from  the  imagery  and  from 
Figure  7. 4.2-1  that  the  two  sets  of  data  had  different  sets  of  sensor  gain  settings,  in  addition  to  other 
envirotunental  differences.  It  was  assumed  that  the  materials  were  relatively  unchanged,  however.  In 
Figure  7.4.2-2  is  also  shown  the  ratio  of  the  post-  to  the  pre-event  “road"  pixels,  which  divided  out  that 
which  was  common  to  the  two  scenes.  The  resultant  pixel  ratio  (ratio  spectrum)  provided  a  correction 
factor  to  compensate  for  differences  in  sensor  gain  factors,  sun  angle  efiects.  atmospheric  path 
attenuation  effects,  and  diffuse  hemispherical  sky  radiance  effects.  Spectral  differences  between  the  road 
and  spoil  are  also  included,  but  they  were  assumed  to  be  minimal  relative  to  the  other  faaors. 

Pixels  were  next  selected  from  spoil  pile  water  ponds  in  the  pre-  and  post-event  scenes. 
Averages  of  nine  pixel  grids  were  generated  from  as  close  as  possible  to  the  identical  water  pond 
locations  in  the  two  images.  The  averaged  pixel  spectra  and  the  ratio  of  the  two  spectra  are  shown  in 
Figure  7. 4.2-3. 
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Figure  7.4.2«1.  Uncorrected  MTL  Before  and  After  Image  Digital  Number  (DN)  Ranges 


FIgur*  7.4.2-2 
RoMl  Pr»*  and  Pott-  Evtnl 


Figur*  7.4^>3 
Wat«r  Pr«>  and  Poat*  Evant 
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Before  dealing  further  with  the  Figure  lA.2-i  spectra  of  the  spoil  pile  water  pond,  additional  .  J 

spectra  are  presented  in  Figures  1A.2-A  and  7.4.2-S.  In  Figure  7.4.2-4  are  average  spectra  of  pixels  frrxn  an 
area  of  the  spoil  pile,  relatively  homogeneous,  away  from  the  ground  zero  detortation  site  and  away 
from  the  cables  and  vehicles.  In  Figure  7.4.2-5  are  average  spectra  of  a  bright  spot,  assumed  to  be  a 
concentration  of  the  fine-grained  fraction,  at  or  close  to  the  ground  zero.  This  spot  was  evident  in  the 
post-event  image  but  not  as  evident  in  the  pre-event  image. 

To  remove  the  spectral  structures  associated  with  the  sensor  gain  factors  and  environmental 
faaors,  the  three  ratio  spectra  (water  pond,  general  spoil,  and  bright  spoO  shown  in  the  bottoms  of 
Figures  7.4.2-3,  7.4.2-4,  and  7.4.2-5  were  divided  by  the  ratio  spectrum  in  Figure  lA.2-2  (the  road  ratio 
spectrum  which  was  used  here  as  a  standard^.  The  resultant  divided  spectra  are  shown  in  Figure  7.4.2-6. 

Examination  of  the  spectra  revealed  unrealistic  discontinuities  in  the  spectra  near  0.44  pm  and  0.84  pm. 

1  he  discontinuities  suggested  the  data  may  have  been  bad  outside  of  the  0.46  to  0.80  pm  wavelength 
range.  Truncated  versions  of  the  spectra  in  Figure  7.4.2-6  are  shown  in  Figure  lA.2-1.  The  speara  in 
Figure  7.4.2-7  represent  the  post-/pre-event  pixel  ratios  for  the  spoil  pile  pond  water,  the  bright  spot,  aitd 
the  spoil.  Dividing  the  original  (water,  bright,  spoil)  post-/pre-event  spectral  by  the  original  road  post- 
/pre-event  (the  standard)  removed  the  gain  factors  and  environment  factors  which  were  common  to 
both  specua  and  produced  the  spectra  in  Figure  lA.2-1,  which  more  accurately  represented  the  "true" 
water,  bright  and  spoil  post-/pre-event  spectra.  ; 

An  additional  ratio  spectrum  is  shown  in  Figure  7.4.2-8  which  represented  the  ratio  of  bright 
spot  material  with  any  residual  road  contribution  removed.  Note  both  the  spoil  pile  pond  water  ratio 
spectrum  (Figure  1A.2~1,  top)  and  the  bright  spot  ra.  general  spoil  ratio  spectrum  (Figure  7.4.Z-8) 
revealed  similar  spectral  ‘reddening*  to  that  observed  in  the  Landsat  TM  data  (Section  7.2). 

The  observation  that  the  reddening  occurred  in  both  the  ASIS  and  TM  imagery,  and  that  the 
pond  reddening  was  stronger  than  the  open  soil  reddening  in  both  datasets  further  supported  the 
possibility  that  the  observed  reddening  was  a  substantive  phenomenon. 
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Figura  7.4.2-7.  Water.  Bright.  Spoi^Road  Pre  and  Post  Event 
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Another  relevant  finding  was  there  were  minimal  fine  scale  spectra  structures  in  the  ASIS  data, 
and  the  large  number  of  spectral  bands  did  not  add  apparent  value  compared  to  the  six-band  Landsat 
TM  data  for  this  observable.  While  high  spectral  resolutions  may  add  value  for  many  applications,  there 
was  no  apparent  advantage  here.  For  this  particular  study  site,  and  this  particular  observable,  the 
spectral  resolution  of  Landsat  TM  revealed  approximately  the  same  level  of  information  as  the  ASIS 
data. 

7.5  THREE-DIMENSIONAL  DIGITAL  MODELING 

For  remote  sensing  projects,  the  use  of  three-dimensional  digital  imagery  models  has  become 
increasingly  common.  These  models  assist  analysts  by  providing  ’three-dimensional’  portrayals  of  a 
given  study  area;  these  models,  however,  are  displayed  in  two-dimensional  image  space.  As  a  result, 
subde  terrain  and  thematic  features  are  shown  in  their  X  and  Y  and  relative  Z  (elevation)  positions. 

The  first  step  in  creating  a  three-dimensional  digital  image  is  to  generate  a  digital  elevation 
model.  Usually,  this  is  accomplished  using  digital  stereo  imagery.  Once  these  elevation  (Z)  data  have 
been  generated,  two-dimensional  (X  and  Y)  imagery  is  registered  to  the  elevation  file.  When  displayed, 
the  incorporated  results  appear  as  a  three-dimensional  image. 

To  generate  the  elevation  dau  of  BEXAR,  stereo  10m  SPOT  panchromatic  images  were 
processed  using  the  Autometric-developed  PEGASUS™  digital  photogrammetric  processing  worksution. 

To  utilize  the  inherent  spatial  resolution  of  the  SPOT  data,  elevation  models  were  generated  at  three 
posting  resolutions.  These  were  esublished  at  10,  20,  and  30m,  respectively. 

An  evaluation  of  the  resultant  models  suggested  the  30m  model,  as  expected,  was  too  coarse.  In 
turn,  the  NTS  elevation  features  were  smoothed  significantly,  and  the  mesa  containing  the  BEXAR  site 
was  not  discernible.  Compared  to  the  30m  generation,  the  20m  model  provided  subsuntially  bener 
definition  of  the  mesa,  and  generally,  the  terrain  shapes  were  distinguishable.  Unfortunately,  the  20m 
generation  represented  a  resolution  threshold  of  the  modeling  process.  Specifically,  the  10m  elevation 
data  did  not  have  an  accepuble  elevation  error  component.  These  findings  were  significant  because 
they  indicated  the  10m  (X  and  Y)  digital  data  could  not  produce  10m  elevation  data.  Rather,  the 
elevation  posting  resolutions  were  dependent  upon  other  modeling  factors  such  as  base-to-height  ratios 
and  relative  off-nadir  locations  of  BEXAR  in  the  two  SPOT  scenes. 

The  next  step  in  the  modeling  process  involved  registering  a  two-dimensional  image  to  the 
elevation  dau.  Both  Landsat  TM  and  SPOT  panchromatic  images  were  accessed  for  this  purpose.  Upon 
displaying  the  resulunt  three-dimensioiul  generation,  however,  it  was  apparent  the  ^m  elevation 
model  was  not  sufficient  for  analyzing  the  BEXAR  event  In  particular,  although  the  mesa  was 
represented  well  in  the  elevation  data,  the  BEXAR  site  —  spatially  comprising  a  portion  of  the  mesa  — 
lost  significant  image  deuil.  Thus,  s^iile  the  general  temin  of  NTS  was  shown  well  using  the  three-  'y 
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dimensional  image  model,  BEXAR,  being  a  site-specific  feature,  was  not  portrayed  sufficiently.  Because 
of  these  results,  no  other  three-dimensional  analyses  were  conducted  using  the  SPOT  elevation  model. 

In  short,  three-dimensional  imagery  models  have  proven  their  utility  in  past  remote  sensing 
research  activities.  Accordingly,  while  the  20m  elevation  model  did  not  produce  adequate  detail  of 
BEXAR,  a  model  with  smaller  posting  resolutions  would  likely  support  site-spedHc  imagery 
investigations  such  as  the  detection  and  monitoring  of  underground  detoiutions. 

7.6  THERMAL  ANALYSES 

Thermal  analyses  of  the  BEXAR  event  were  conduaed  using  Laixlsat  TM  Band  6  imagery.  In 
general.  Band  6  (10.4  to  12.5  pm)  measures  the  amount  of  infrared  radiant  flux  emitted  from  surfaces. 
Apparent  temperature  is  a  function  of  the  emissivity  and  true  kinetic  temperature  of  a  given  surface. 
Initially,  the  thermal  dau,  collected  before  and  after  the  event,  were  to  be  analyzed  for  detections  of 
event-related  change.  Upon  displaying  the  thermal  imagery,  however,  it  was  unfortunately  apparent  the 
120m  spatial  resolution  of  Band  6  was  not  adequate  for  interpreting  change.  In  short,  while  gerteralized 
NTS  seasonal  thermal  changes  were  observed,  because  of  the  p<x)r  spatial  resolution,  no  site-specific 
thermal  change  related  to  Che  BEXAR  evetu  was  found 

Surface  thermal  change  remains  a  potential  observable  for  the  detection  and  monitoring  of 
underground  detonations.  For  such  change  to  be  observed,  though,  the  thermal  imagery  must  have 
significandy  higher  spatial  resolutions  than  that  provided  by  the  120m  TM  Band  6  data.  Likely,  spatial 
resolutions  in  Che  10m  range  would  be  reciuired 
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8.1 


INTRODUCTION 


SECnON  8 

INTKODUCnON/BACKGROUND 


The  purpose  of  this  portion  of  the  research  was  to  define  an  integrated  workstation  concept 
addressing  the  requirements  for  economical,  timely,  reliable,  and  versatile  exploitation.  Indeed,  the 
increasing  capabilities  of  medium-  to  low-cost  workstations  has  allowed  'commercal-off-the-shelf ,  or 
COTS,  technologies  and  standards  to  tie  applied  to  newly-emerging  analytical  areas.  A  system  was 
required  to  support  end-to-end  solutions  (i.e.,  data  import,  processing,  and  output)  and  be  designed 
based  on  industry  standards  (e.g.,  UNIX,  C,  W-Windows).  The  combination  of  UNIX,  C,  W-Windows, 
standardized  bus  architectures,  and  general  purpose  processors  (GP)  provides  an  analyst  with  a 
processing  platform  commonly  referred  to  as  an  cpen  system.  The  open  system  concept  is  discussed 
here  for  its  role  in  detecting,  monitoring,  and  analyzing  urKierground  nuclear  tests  and  associated  test 
facilities. 

The  consolidation  and  stabilization  of  operating  systems,  software  development  languages,  and 
graphical  user  interfaces  has  allowed  software  development  activities  to  focus  on  the  development  of 
more  modular  and  sophisticated  software  packages.  Database  generation  and  maintenance  aaivity  is 
required  for  many  applications,  and  because  of  this,  it  can  be  considered  an  enabling  technology  (those 
portions  of  software  that  application  developers  use  as  part  of  their  total  solution).  In  contrast, 
applications  provide  a  total  solution  to  a  user's  problem,  such  as  nuclear  test  monitoring,  targeting, 
indications  and  warning,  mission  rehearsal,  cross-country  mobility  analysis,  etc. 

8.2  SYSTEM  FUNCTIONING 


The  fusion  of  multiple  products  to  produce  a  working  display  for  detecting  and  monitoring  test 
facilities  will  demand  unique  software  and  dedicated  equipment  Fundamental  to  solvirtg  the  data  fusion 
problem  will  be  the  ability  to  bring  multisensor  imagery  having  different  viewing  geometries  into  a 
common  and  accurate  coordinate  system.  Disparate  datasets  can  be  registered  to  each  other  directly  or 
indirecdy  by  being  first  tied  to  geo-coordinates  of  the  earth.  Once  a  merged  dau  set  has  been 
generated,  many  analysis  options  become  possible.  For  nuclear  test  detection  and  monitoring  scenarios, 
three  applications  were  considered  as  important  stereo  viewing  and  mensuration,  three-dimensional 
terrain  visualization  and  analysis,  and  multispectral/hyperspectral  analysis.  Spedflcally,  the  exploitation 
concept  provided  includes  the  following  capabilities; 


) 
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Digitizing  hardcopy  daU; 

Abilities  to  fuse  various  data  types; 

Base  imagj  enhancement; 

Multispectral/Hyperspectral  analysis; 

Precise  Mensuration; 

Stereo  viewing; 

Three-dimensional  perspective  scene  generation;  and 

High-quality  hardcopy  output. 

Figure  8.2-1  illustrates  the  major  functional  components  of  a  digital  exploitation  workstation. 
There  are  two  basic  components  to  the  architecture.  The  first  comporrent  (shown  on  the  left  portion  of 
Figure  8.2-1)  is  Database  Generation  and  Maintenance;  this  includes:  (1)  the  importation  of  digital  and 
hardcopy  datasets  (input  functions);  and  (2)  generation  of  merged  datasets  using  a  variety  of 
photogrammetric  techniques.  Rectified  imagery  can  be  used  to  generate  enhancements  to  the  database 
such  as  feature  data,  digital  elevation  models,  or  orthorectified  images.  These  data  can  be  used  for 
applications  (second  basic  component)  shown  on  the  right  of  Figure  8.2-1.  By  keeping  the  data 
preparation  software  separate  from  the  applicatiorts  software,  duplication  of  an  overlap  between  software 
functions  is  reduced.  This  also  enables  rapid  prototyping  of  new  applications  oriented  to  specific 
analytical  requirements.  The  various  system  components  are  discussed  separately  in  the  following 
sections. 

8.3  DATABASE  GENERATION  AND  MAINTENANCE 

8.3.1  Hardcopy  Data  Import 

Hardcopy  dau  must  be  converted  to  a  digital  format  to  be  utilized  by  a  digital  exploitation 
system.  This  usually  is  accomplished  with  a  digitizer  (also  known  as  a  scanner).  The  type  of  digitizer 
required  for  hardcopy  input  is  defined,  in  large  part,  by  the  type  of  hardcopy  dau  being  imported  to 
the  image  processing  system.  If  map  daU  are  to  be  incorporated,  an  inexpensive  flatbed  scanner  with  a 
resolution  of  300  dots/inch  (dpi)  can  be  used  (approximate  cost  is  $7,0(X)).  If  imagery  is  to  be  digitized, 
a  higher  resolution  system  is  usually  required  such  as  an  Eikonix  1412  digitizing  camera  with  a  Gordon 
Instmments  digitizing  worksution.  These  can  digitize  4000-by-40(X)-pixel  scenes  from  black  and  white  or 
color  prints  and  transparencies  at  a  resolution  of  up  to  20  micrometers  (approximate  cost  is  $30,0(X)).  .  > 
high-end  option  exists  where  a  high-resolution  digitizing  capability  is  requited  over  large  areas.  In  this 
scenario,  a  digitizer  with  a  movable  suge  is  used  to  move  accurately  around  large  images  arid  digitize 
smaller  patches.  These  are  then  mosaicked  together  (approximate  cost  is  $60,000). 
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8.3.2  Softcopy  Dau  Import 


There  are  an  increasing  number  of  ways  to  import  digital  data  to  digital  exploitation  systems. 
Most  digiul  imagery  is  stored  on  magnetic  tape  media.  Tape  formats  include:  150  Mbyte  SCSI  cassettes, 
common  on  many  engineering  workstations;  9-track  tapes  which  hold  approximately  the  same  amount 
of  dau  as  SCSI  cassettes;  and  a  recently-developed  arxi  inexpensive  upe  storage  capability  using  8mm 
upes  (SGbytes). 

Another  recently  introduced  format  and  media,  which  are  becoming  increasingly  popular  to 
users,  are  CD-ROMs  produced  by  the  Defense  Mapping  Agency  (DMA).  The  disks  are  similar  to  CDs 
commonly  used  in  home  stereos,  and  these  can  hold  large  amounts  of  dau  (650  Mbytes.).  DMA 
currently  distributes  a  number  of  dau  types  on  CD-ROMs  including  Arc  Digiul  Raster  Graphic  (ADRG) 
dau  (digiul  maps).  Digital  Terrain  Elevation  Dau  (DTED),  and  some  vector  datasets  such  as  World 
Vector  Shoreline  (WVS).  The  capability  to  read  these  new  datasets  is  critical  to  an  imagery  exploiution 
system. 

8.3  3  Data  Storage 

The  various  dau  utilized  by  digital  exploiution  systems  require  large  amounts  of  dau  storage. 
Dau  storage  remains  a  significant  problem  for  digital  exploiution  systems.  A  standard  three-color  9-  by 
9-inch  aerial  image,  digitized  at  100  dpi  resolution,  conuins  approximately  2.45  *  10^  pixels.  Most 
remote  sensing  projects  contain  many  such  large  images  which  must  be  stored  in  a  system  during  the 
compilation  process.  If  they  can  be  read  into  the  system  from  magnetic  upe,  the  upe  reading  itself 
takes  considerable  time. 

Generally  hardware  vendors  have  not  yet  provided  dau  storage  devices  to  hold  several 
gigabytes  of  information  cost-effectively.  Tape  cartridge  storage  is  inexpensive  but  requires  significant 
seek  and  transfer  time  (more  than  ten  minutes  to  locate  one-gigabyte  images).  While  dau  storage  costs 
are  continuously  dropping,  daU  storage  will  remain  a  major  bottleneck  in  projects  using  large  volumes 
of  softcopy  imagery. 

With  regard  to  format  coiuiderations  as  much  as  possible,  common  dau  formau  should  be 
used.  For  insunce,  the  dau  file  format  may  be  a  generic  'flat  Hie*  with  neither  header  information  nor 
proprieury  formatting.  Alternatively,  it  may  be  in  a  proprietary  format  of  some  commercial  relational 
daubase  management  system  (RDBMS).  Other  proprieury  or  special  formau  may  be  required  to 
support  COTS  or  custom  application  software  packages. 

If  multiple  formau  exist,  to  maruge  the  dau  Hies  shared  between  the  various  software 
applicatiotu  and  the  woriutations  die  ideal  solution  would  be  a  virtual  database  inietftne  (VIM).  The  VDI 
serves  as  the  clearinghouse  for  all  dau  exchai^.  Autometric  is  striving  toward  developirtg  such  a 
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solution;  however,  it  is  not  a  COTS  solution  at  this  time.  If  a  VDI  is  not  available,  software  utilities  should 
be  generated  to  convert  between  the  various  data  formats  required  by  the  many  applications. 

8.4  APPLICATIONS 

8.4.1  Data  Rectification 

Without  a  method  for  registration  of  imagery  to  a  geographic  coordinate  system  or  other 
images  in  a  rigorous,  accurate  manner,  mensuration  capabilities  are  in  doubt.  Often  vector  overlay  dau 
and  elevation  dau  are  tied  to  the  earth;  and  without  knowledge  of  where  images  are  located  on  the 
earth,  it  is  difficult,  if  not  impossible,  to  fuse  the  datasets. 

Over  the  last  30  years.  Autometric  has  based  a  significant  portion  of  its  business  on  its  often 
unique  capability  to  relate  imagery  to  accurate  geo-positions.  The  heart  of  this  capability  is  the  use  of 
rigorous  and  accurate  imaging  sensor  models  that  provide  for  simuluneous  registrations  of  up  to  ten 
overlapping  or  independent  images  to  ground  control.  Autometric  has  emphasized  accuracy  of 
mathematical  models  and  dau  integrity  in  the  development  of  its  dau  fusion  and  mensuration 
capabilities  (particularly  since  much  of  Autometric's  work  has  involved  defense  mapping,  miliury 
targeting,  and  intelligence  exploiution).  Miliury  sensors  frequently  have  highly  complex  imaging 
characteristics  which  require  rigorous  adherence  to  the  modeling  of  the  physical  phenomena  of  the 
imaging  event  to  ensure  the  mathematical  models  correctly  represent  the  sensor.  The  mathematical 
models  present  in  the  Autometric-developed  photo^mmetric  software  preserve  the  metric  accuracies 
inherent  in  these  complex  imaging  systems. 

Once  a  reaifled  imagery  daubase  has  been  constructed,  the  mensuration  functionality  is  not 
complicated.  Using  polynomials  generated  by  the  sensor  models,  geographic  coordinates  can  be 
determined  for  any  point  in  the  image  allowing  disunces,  azimuths,  and  areas  to  be  calculated  simply 
and  accurately. 

8.4.2  Multispectral  and  Hyperspectral  Image  Processing 

The  most  important  hardware  consideration  when  processing  multispectral  dau  is  the  display 
device;  24-bit  color  is  an  option  on  most  systems  but  not  necessarily  a  standard  configuration. 
Additiorully,  if  imagery  is  to  be  scrolled  or  three-dimensional  inuges  are  to  be  manipulated,  a  double 
buffering  capability  may  be  necessary,  and  this  requites  48-bits  of  ct^or. 

As  for  tire  software,  basic  multispectral  analysis  capabilities  ate  readily  available  on  a  wide 
variety  of  hardware  conf^ratiotts.  It  is  important,  however,  to  determine  tint  die  softarare  package 
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selected  can  accept  and  process  the  expected  input  dau.  For  example,  hyperspectral  analysis 
capabilities  are  not  as  commoa 

Besides  COTS  solutions,  there  are  also  custom  processing  techniques  available  which  hold  great 
promise  for  analyses  of  multispecual  and  hyperspectral  imagery  (e.g.,  subpixel  analyses). 

Given  the  map-like  quality  of  commercial  multispectral  imagery,  users  often  perform  analyses 
on  these  dau  using  other  non-imagery  dau  such  as  elevation  dau  and  thematic  (map)  dau.  These 
capabilities  are  provided  in  geographic  information  systems  (GIS).  Conuiiung  daU  maiugement 
functions,  these  are  slightly  more  complex  than  standard  multispectral  image  processing  packages.  It  is 
impoitant  that  the  image  processing  and  GIS  sctftware  (usually  separate  packages)  utilized  allows  for 
exchange  of  dau  between  these  two  functional  areas. 

8.4.3  Stereo  Viewing  and  Mensuration 

With  ceruin  sensors  (e.g.,  SPOT),  stereo  pair  images  can  be  acquired.  Using  stereo  viewing 
techniques,  the  terrain  is  depicted  in  three  dimensions  which  is  extremely  useful  for  terrain, 
topographic,  and  microchromatological  arulyses.  Prior  to  viewing,  the  stereo  pairs  must  be  prepared 
using  rectification  procedures. 

In  addition  to  viewing  in  stereo,  a  complimenury  capability  is  the  ability  to  mensurate 
accurately  in  three  dimensions.  If  the  imagery  has  been  tied  to  the  earth's  surface,  the  latitude,  longitude, 
and  elevation  (or  similar  coordinates)  can  be  determined  for  any  point  in  the  image.  This  is  an 
important  capability  when  measuring  locations,  aerial  extent,  and  volumes  of  various  features  found  in 
the  imagery. 

8.4.4  Three-Dimensional  Terrain  Visualization 

Stereo  imagery  nuy  ikx  always  be  available  for  an  area  of  interest  Alternatively,  the  terrain  can 
still  be  modeled  in  three  dimensions  if  imagery  and  elevation  daU  (e.g.,  DTED  dau)  ate  both  available. 
Using  perspective  transformation  techniques,  an  accurate  represenution  of  the  terrain  can  be  created. 
On  some  hardware  platforms,  sudi  as  most  Silicon  Graphics  worksUtions,  perspective  scenes  can  be 
built  quickly  and  can  be  manipulated  interactively. 

8.4.5  Applications  Developmem  Environment 

In  addition  to  using  COTS  software  paduges,  a  digital  exploitation  wortcstation  should  provide 
an  applicatiom  developmertt  environment.  Industry  sundards  such  X-Windowa,  UNIX,  MOTIF,  arxl  C 
should  be  supported.  When  available,  development  optioiu  should  be  ao^iied  fur  the  COTS 
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applications.  These  options  usually  consist  of  libraries  of  software  funciions  that  can  be  used  to  add 
functionality  to  the  COTS  software  or  to  build  separate  applications. 

8.5  HARDCOPY  OUTPUT 

The  majority  of  analyses  of  digital  imagery  will  be  accomplished  in  softcopy;  however,  there  is 
always  a  need  for  the  production  of  hardcopy  output  Hardcopy  output  is  necessary  to  record  images  as 
displayed  on  a  high-resolution  computer  display.  Hardcopy  devices  can  produce  output  that  contains  a 
larger  image  area  than  can  be  displayed  on  a  computer  display  screen.  Printed  output  can  be  used  as  a 
storage  medium  and  can  be  easily  utilized  by  many  people  without  specialized  computer  equipment  For 
this  reason,  and  others,  results  from  softcopy  analyses  are  often  presented  in  a  hardcopy  format 

V(^n  working  with  multispectral  image  dau,  it  is  important  that  the  printing  device  be  capable 
of  24-bit  "true-colof*  printing.  In  many  cases,  three  image  spectral  bands  are  combined  and  printed  as  a 
red-green-blue  (RGB)  image.  Color  printers  are  capable  of  printing  black  and  white  images  by  repeating 
the  same  data  to  the  three  RGB  charmels  of  the  printer.  However,  sometimes  this  method  produces 
black  and  white  images  that  have  a  color  cast  to  them  (e.g.,  bluish  tinge).  If  only  black  atKl  white  images 
are  to  be  produced,  then  a  black  and  white  printer  should  be  used  or  the  color  printer  should  be 
adjusted  accordingly.  Depending  on  the  type  of  printing  process,  this  adjustment  might  include  using  a 
black  and  white  ribbon  (thermal  dye  printer)  or  black  ink  exclusively  (ink-jet  printer). 

There  are  numerous  hardcopy  output  options  available.  The  one  seleaed  depends  on  the 
hardcopy  output  needs.  If  high  quality  output  is  desired  (photographic  quality),  a  continuous-tone 
process  should  be  used.  The  3M  Color  Laser  Imager  (CLI)  utilizes  a  wet  photographic  process  that 
produces  extremely  high-quality  color  output  at  300  pixels/inch  resolution.  Its  output  size  is  limited  to  12 
by  18  inches.  This  printer  is  a  high-throughput  model  that  costs  approximately  $133,000.  For  lower 
quantity  and  smaller  format  printing,  the  Kodak  XL7700  printer  is  a  good  alternative.  Its  image  quality  is 
not  as  high  as  the  CU  since  it  records  at  200  pixels/inch.  It  utilizes  a  thermal  dye  process,  and  its 
throughput  is  not  as  high  as  the  CLL  The  XL7700  output  is  limited  to  11  by  11  inches.  The  ct>st  of  a  color 
XL7700  is  approximately  $20,000. 

If  larger  format  output  is  desired,  contiimous  tone  printing  is  not  viable;  however,  there  are  three 
options  available.  The  first  option  is  photographic  enlargement.  This  process  is  expensive,  and  if 
geometric  accuracy  of  the  data  is  to  be  preserved,  this  is  ttot  a  realistic  option.  Litht^raphic  printing 
can  be  done,  but  its  cost  is  only  justified  vdien  large  numbers  of  prints  are  to  be  produced  (minimum  of 
500).  The  lithographic  process  utilizes  much  specialized  equipment,  and  the  image  quality  is  inferior  to 
continuous  tone  printing  since  half-tones  are  used.  The  cost  of  a  medium-fonnat  UtlK>graphic  image 
setter  arxl  pre-press  erjuipment  would  be  approximately  $130,000,  and  this  would  not  indude  printing 
press  plate  production  equipment.  The  third  large  format  option  is  die  DUS  Graphics  3024  or  3047 


printer.  These  two  printers  are  ink-fet  printers  that  offer  exceptional  image  quality  at  300  pixels/inch 
resolution.  To  the  unaided  eye,  the  output  appears  to  be  a  continuous  tone;  however,  these  printers 
actually  utilize  fine  resolution  half-tones.  For  most  purposes,  these  half-tones  are  not  a  hindrance,  unless 
the  image  data  are  to  be  aiulyzed  under  magnification.  The  output  sizes  for  these  two  printers  is  24  by 
24  inches  for  the  3024  and  34  by  47  itKhes  for  the  3047.  Therefore,  high-quality  full-size  image  maps  can 
be  produced  on  the  3047.  The  3047  printer  sells  for  $130,000  ar>d  the  3024  for  $84,500. 

8.6  SYSTEM  CONCEPT 


The  rationale  used  for  the  analysis  system  is  to  develop  an  open  system  built  using  COTS 
solutions  as  much  as  possible:  this  includes  both  hardware  and  software.  The  major  advantage  to  this 
approach  is  that  it  gready  reduces  costs  (both  initial  and  maintenance).  A  typical  digital  exploitation 
worksution  is  shown  in  Figure  8.6-1. 

The  suggested  host  computer  is  a  Silicon  Graphics  4D/35TG  workstation.  The  UNIX-based 
workstation  has  a  COTS  general  purpose  36  Mhz  RISC  cpu.  It  utilizes  COTS  graphics  boards  that  enable 
24-bit  color  display  of  imagery  and  also  double-buffering  of  12-bit  data.  It  easily  handles  large  image 
dausets  quickly.  The  worksution  can  also  support  a  hardcopy  digitizer  and  hardcopy  printer.  The 
Silicon  Graphics  is  a  COTS  hardware  solution  that  supports  industry  standards  for  software  development 
with  UNIX,  X-Wirxlows,  MOTIF,  and  C,  to  name  a  few. 
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Since  the  Silicon  Graphics  computer  is  a  COTS  solution,  there  is  a  variety  of  COTS  software 
available  to  support  this  platform.  There  are  various  options  to  consider  for  both  image  processing,  GIS, 
terrain  visualization,  and  photogrammetiic  (rectification  and  mensuration)  functionality.  As  an  analysis 
laboratory  grows  in  size,  additional  workstations  can  be  added  via  a  Local  Area  Network  (LAN),  and  the 
dau  files  can  shared. 

For  hardcopy  digitization,  the  Eikonix  1412  scarmer  is  suggested.  This  scarmer  along  with  the 
Gordon  Instruments  transmissive/reflective  workstation  will  provide  capability  to  digitize  black  and 
white  and  color  paper  and  transparent  imagery. 

The  most  viable  solution  for  small-format,  high-quality  printing  is  the  Kodak  XL7700  printer. 
This  printer  is  easy  to  use  and  can  make  both  black  and  white  and  color  prints  (separate  ribbons).  An 
additional  printer  may  be  desired  depending  on  the  required  type  of  output  A  decision  should  be 
made  between  a  high-throughput  CU  that  produces  12  by  18  inch  output  or  a  large  format  printer  such 
as  the  IRIS  Graphics  3047  (which  may  take  an  hour  to  produce  a  34  inch  by  47  inch  color  prinO- 

8.7  COSTS 

Costs  are  provided  for  workstation  components  discussed  in  Section  8.6.  These  are  approximate 
costs;  no  discounts  have  been  considered. 

Silicon  Graphics  4D/35TG  Workstations  $54,000 

48  Mb  RAM 
1.2  Gbyte  hard  disk 
9-irack  tape  drive 
8  mm  cassette  drive 
CD-ROM  drive 

Eikonix  1412  Digitizer  and  Gordon 


Instruments  Digitizing  Worksution 

Application  Software  Packages 

Qnchiding  GIS,  image  processing, 

$30,000 

photogrammetry,  terrain  visualization) 

$90,000 

Kodak  XL7700  Color  Printer 

$20,000 

Printer  Options: 

IRIS  Graphics  3024  Printer 

$84,500 

IRIS  Graphics  3047  Primer 

$130,000 

3M  Ctrfor  Laser  Imager 

$133,000 
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SEcmcm  9 
CONCLUSIONS 

9.1  STUDY  CONCLUSIONS 

Based  upon  the  analyses  conducted  under  this  two  phase  research  study,  the  following 

conclusions  are  documented: 

9.1.1  General 

•  With  the  dissolution  of  the  Soviet  Union  and  proliferation  of  nuclear  technology,  the 
requirement  for  global  search,  detection  and  monitoring  of  current  and  pending  nuclear  test 
facilities  increases  in  priority. 

•  The  use  of  current  and  planned  national  and  international  civil  and  commercial  imaging 
satellites,  combined  with  National  assets,  provide  a  reliable  means  for  detecting  and  monitoring 
worldwide  underground  nuclear  test  preparations  and  detonations,  whether  covert  or  under 
arms  control  treaties. 

•  The  large  area  coverage  and  wide  area  search  capabilities  of  multispectral  satellite  systems  such 
as  Landsat  and  SPOT  would  aid  in  establishing  whether  nuclear  weapons  development  has 
advanced  beyond  the  laboratory,  by  providing  important  information  on  test  site  activities. 

•  The  stability  of  spacebome  platforms  for  imaging  systems  is  far  superior  to  that  of  aircraft 
platforms  which,  and  when  combined  with  their  synoptic,  repetitive,  wide  area  coverage 
parameters,  results  in  the  conclusion  that  spacebome  platforms  enhance  exploitation  and,  in 
particular,  the  registration  of  scenes  for  change  detection. 

9.1.2  Database  Management 

•  Establishing  a  database  of  remote  sensor/reconnaissance  imagery,  both  hard  and  softcopy,  is  a 
key  componeiu  for  monitoring  worldwide  underground  nuclear  develo{»nent  and  testing. 

•  Once  an  area  becomes  suspect,  all  source  historical  data  of  that  area  should  be  gathered  to 
esublish  a  comprehensive  daubase  that  will  support  collection  planning  and  product 
exploitation. 

•  The  collation  of  both  imaging  sensor  products  and  coUateial  dau  such  as  maps  and  schematics 
into  a  common  digital  media  is  the  most  efficient  means  for  exptoidng  the  multiple  dau  sets. 
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•  Generic  test  preparation  chronological  timelines  can  be  documented  through  the  use  of 
commercial  imaging  satellites,  and  can  provide  a  reliable  means  for  monitoring  the  progress 
and  eventual  success  of  future  underground  tests. 

9.1.3  Collection  Management 

•  "Smart”  all-source  imagery  collection  planning  is  mandatory,  both  in  terms  of  global  and 
directed  search  strategies  and  in  selective  multiple  sensor/multispectral  collections,  once  a 
suspicious  area  is  delineated. 

•  Periodic  and  point  target  tasking  of  imaging  sensors  should  begin  at  the  inception  of  a  test  and 
continue  past  the  day  of  the  final  event. 

•  Collection  management  and  tasking  of  imaging  sensors  should  be  coordinated  with  other 
intelligence  and  verification  sources,  including  seismic  monitoring. 

9.1.4  Data  Processing/Exploiiation 

•  The  combined  use  of  various  imaging  systems  and  the  correlation  of  multiple 
sensor/multispectral  imagery  with  collateral  dau  demands  the  use  of  a  computer  based  digital 
image  processing,  exploiution  and  fusion  system. 

•  A  test  bed  digital  workstation  designed  and  dedicated  to  the  requirements  for  detecting  and 
monitoring  underground  tests  should  be  developed  and  installed  within  the  DARPA  and 
eventually  deployed  operationally. 

•  The  30m  spatial  resolution  of  the  Landsaf  TM  is  marginal  in  terms  of  detecting  activity  at  a  site. 
The  15m  spatial  resolution  of  the  pending  Landsat  6  panchromatic  band  will  allow  the  detection 
of  heavy  equipment  and  trailers  at  a  site. 

•  Ground  scarring  and  terrain  changes  can  best  be  analyzed  using  the  Landsat  TM  Band  5  and  7 
spectral  windows  whereas  vegeutive  changes  are  emphasized  by  using  the  visible  and  near 
infrared  bands. 

•  Tools  have  been  developed  and  tested  in  this  study  that  will  allow  the  isolation  of  changes  that 
occur  due  to  testing  from  those  that  result  from  nominal  atmospheric  or  vegeutive  changes. 

•  A  precise  digiul  thiee-dimensional  perspective  view  of  an  identified  test  area  should  be 
produced  to  support  analyses  of  the  site,  the  future  tasking  of  collection  asseu,  and  a  readily 
understood  vehicle  for  dedsioiunakers. 

9.2  BEXAR  TEST  CONCLUSIONS 
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9.2.1  General 


•  For  a  given  test  having  a  spoil  pile  at  or  near  ground-zero,  the  detonation  should  change  the 
pile  predictably.  Phenomenologically,  the  detonation  event  should  produce  a  variety  of  observables. 
From  the  research,  these  changes  were  categorized  as  either  spectral  changes  (e.g.,  reddening,  cracking 
and  shadowing,  particulate  shifts)  and/or  spatial  changes  (e.g.,  lateral  expansion).  The  usefulness  of 
these  indicators,  however,  would  likely  be  dependent  on  the  size  and  depth  of  the  spoil  pile,  depth  of 
the  detonation,  and  the  characteristics  of  the  overburden. 

9.2.2  Multispectral 

•  BEXAR-related  changes  were  found  using  the  multispectral  TM  data,  and  these  were  detected 
using  statistical  techniques  to  model  the  spectral  variance  measures.  For  detonation  conditions 
comparable  to  those  of  BEXAR,  the  results  suggested  Landsat  TM,  at  30m,  provided  adequate  spatial 
resolution  capabilities  for  detecting  event-related  changes.  Furthermore,  the  statistical  techniques 
required  minimal  complementary  visual  analyses  of  the  data. 

•  Soil  and  geologic  changes  detected  were  attributed  to  the  BEXAR  event 

•  Vegetation  changes  were  due  to  naturally-occurring  plant  cydical  states  and  were  not  attributed 
to  the  BEXAR  event 

•  Originally,  the  TM  spectral  band  positions  and  the  bandwidths  were  designed  to  be  effective  for 
a  multitude  of  diverse  earth  resource  monitoring  tasks.  These  spectral  bands  were  sufficient  for 
underground  detonation  detections  because  they  addressed  the  diverse  observables  associated  with  the 
test 

•  A  space  borne  platform,  such  as  Landsat,  is  best  for  detecting  change  because  of  its  orbiul 
stability.  Images  can  be  acquired  repeatedly  of  the  same  area  from  approximately  the  same  position  in 
space.  In  particular,  the  subility  of  Landsat  is  superior  to  that  of  an  airborne  platform  subject  to  pitch, 
yaw,  and  roll  distortions. 

9.2.3  Subpixel 

•  BEXAR-related  changes  were  found  using  the  subpixel  technique.  Specifically,  the  process 
revealed  the  occurrence  of  lateral  expansioas  of  the  BEXAR  site's  surface. 

•  Vegetative  rhangt»«  were  deiecKd,  but  none  of  diese  were  attrfixited  to  the  KXAR  event 

•  The  surface  expansion  was  defined  better  using  the  subpixel  technics  compared  to  the  whole- 
pixel  analyst  of  spectral  variatuse  mcthoda;  and  accordingly,  the  subpixd  process  oooifrfhiiented  well 
the  multispeotal  tedmiques. 
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•  Post-event  spectral  reddening  of  the  site-specific  iron  oxide  found  using  the  whole-pixel 
multispectral  analyses  was  also  found  using  the  subpixel  analyses  technique. 

9.2.4  Hyperspectral 

•  BEXAR-related  changes  were  found  using  the  ASIS  hyperspectral  data.  Specifically,  discriminate 
soil/surface  changes  (non-statistical)  were  attributed  to  the  BEXAR  event;  however,  no  event-related 
vegetative  changes  were  found. 

•  Due  to  specific  sensor  malfunctions  (in  the  ASIS  specuometer  number  three),  the  greater 
spectral  resolution  provided  by  ASIS  did  not  provide  substantially  more  information  vs.  the  six-band 
TM  data;  greater  spatial  resolutions,  however,  did  allow  for  detection  of  change  not  found  using  the  TM. 

•  From  an  imagery  perspective,  the  ASIS  data  were  not  used  in  their  original  l6-bit  radiometric 
range,  and  subsequently,  information  was  lost 

•  For  small,  site-specific  applications,  the  utility  of  airborne  sensor  data,  such  as  ASIS,  is  limited 
due  to  significant  impacts  of  aircraft  pitch,  roll,  and  yaw  distortions. 

9.2.5  Thermal 

•  Due  to  the  poor  spatial  resolution  of  the  TM  Band  6  (120  by  120m),  no  thermal  changes 
attributable  to  the  BEXAR  event  were  found. 

9.2.6  Digital  Three-Dimensional  Modeling 

•  Digital  elevation  data  generated  using  10m  SPOT  stereo  imagery  did  not  provide  sufficient 
spatial  definition  clarity  for  small-area  site-specific  modelings  of  the  BEXAR  site. 


10.1  STUDY  RECOMMENDATIONS 


SECTION  10 
RECOMMENDATIONS 


10.1.1  General 

•  An  ideal  sensor  configuration  and  acquisition  strategy  should  be  considered  for  detecting  test- 
site  potential  and  actual  observables.  It  is  recommended  that  observables  at  NTS  and  other  test  sites  be 
evaluated  further  before  specific  sensor  and  sensor  parameters  are  declared.  Indeed,  other  observables 
should  be  verified  as  representative  of  event-related  change  phenomena.  Detectabilities  of  these 
additional  observables  using  LarKlsat  or  other  airborne  sensors  must  be  determined. 

•  Many  of  the  environmenul  changes  due  to  the  BEXAR  blast  were  subtle.  Some  of  the  changes 
were  evident  immediately  after  the  blast  whereas  other  changes  were  seemingly  not  apparent  for  days. 
Other  potential  changes  could  take  weeks  or  months  to  occur  (such  as  vegetation  stress).  For  these 
reasons,  an  image  collection  rate  of  approximately  once  every  seven  days  would  be  desirable  for  a 
given  test  site.  The  seven-day  rate  would  also  help  quantify  natural  rates  of  environmental  change. 
(These  can  be  quite  rapid  depertding  on  the  growth  season  status.) 

•  An  economical  digital  worksution  should  be  developed  that  could  be  tested  and  deployed  to 
field  sites  to  support  treaty  verification  and  site  monitoring. 

10.1.2  Multispectral 

•  Landsat  TM  should  be  considered  as  a  primary  sensor  for  monitoring  and  detecting 
underground  explosions.  This  is  due  to  Landsafs  coUection-to-coUection  geometric  fidelity.  Thus, 
Landsat  is  not  only  capable  of  sealing  event-related  changes,  but  the  TM  dau  could  serve  as  the 
geometrically-accurate  registration  base  for  subsequent  multisensor  fusions  (patticulariy  when  adequate 
map  data  are  not  available). 

•  Research  should  be  focused  on  the  potential  utility  of  added  spatial  resolution  capabilities  of 
future  Landsat's  (i.e.,  6,  7,  and  8) 

•  To  provide  more  quantificationa  of  change,  analysis  of  variance  techniques  should  be  used  in 
conjunction  with  advanced  regression  techniques  (e.g.,  USREL)  to  measure  latent  Cumneasured”) 
variables  of  a  given  site. 
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10.1.3 


Subpixel 


•  Since  the  subpixel  process  can  detect  objects  comprising  as  little  as  20%  of  a  spatial  pixel,  the 
utility  of  the  technique  should  be  researched  in  conjunction  with  the  spatial  resolution  changes  of 
pending  multispectral  systems. 

•  Since  the  subpixel  process  requires  at  least  five  multispectral  bands  and  accuracy  increases 
when  more  bands  are  available,  the  utility  of  the  technique  should  be  researched  in  conjunction  with 
the  spectral  resolution  changes  of  pending  multispectral  systems. 

•  Other  advanced  subpixel  techniques  exist,  and  their  validities  and  reliabilities  should  be 
evaluated  for  purposes  of  detecting  and  monitoriiig  underground  tests. 

•  Due  to  the  computational  demands  of  subpixel  processings,  high-end  computing  devices  — 
such  as  the  NIDL  Princeton  Engine  —  should  be  evaluated  for  their  potential  application  and  testing 
roles. 

10.1.4  Hyperspectral 

•  Hyperspectral  data  must  be  analyzed  both  as  image  data  and  non-image  data;  specifically,  l6-bit 
time-series  and  frequency-domain  stochastic  modeling  techniques  should  be  considered  as 
complements  to  8-bit  discriminate  imagery  processing. 

•  If  aircraft  hyperspectral  dau  are  to  be  used,  collection  management  parameters  must  be  strictly 
controlled. 

•  As  other  hyperspectral  and  ultraspectral  systems  evolve,  their  comparative  utilities  for  detecting 
and  monitoring  underground  tests  must  be  considered. 

10.1.5  Thermal 

•  The  spatial  resolution  of  the  thermal  data  must  be  signihcantly  better  than  the  120m  Landsat  TM 
thermal  imagery.  Likely,  a  spatial  resolution  in  the  10m  range  would  be  needed  if  thermal  data  are  to 
detea  event-related  changes. 

•  As  multiband  thermal  systems  evolve,  their  comparative  utilities  for  detecting  and  monitoring 
underground  tests  must  be  considered 

10. 1.6  Digital  Three-Oimertsional  Modeling 

•  if  digiul  three-dimensional  models  of  a  given  test  site  are  to  be  generated,  the  stereo  dau  must 
have  sigruficantly  better  spatial  resolutions  than  diat  offered  by  the  10m  SPOT  system. 
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•  A  precise  digital  thiee-dimensi^  nal  perspective  view  of  an  identified  test  should  be  provided  to 
support  analyses  of  the  site,  future  tasking  of  collection  assets,  and  decision-maker  and  analyst  viewing. 
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